NASA Technical Memorandum 102467 


( N A S 
J! ->I 
l N A c - 



■1 

1 


FLUSH: A Tool for the Design of Slush 
Hydrogen Flow Systems 




Terry L. Hardy 
Lewis Research Center 
Cleveland, Ohio 


February 1990 


NASA 

A- T M - 1 > 7 ) i-LU'X: A TU^l H)K T Ht N r > 0 -l 7 i >90 

,N ’’f AU' : 1 1 HV-MG-N r L n W 

A) /I L SCI ‘’I! 

Unci -is 

>->/<?'< 0/>oAjLj 


-5255 


^ ****m*"fm* m * ' * 'Ti >rr -r |imn»[i.inilW»JWIf»>IL~IW|-|UI«-inili .. tv^Ti-T-n- n-irr ;- i 


FLUSH: A TOOL FOR THE DESIGN OF SLUSH HYDROGEN FLOW SYSTEMS 

Terry L. Hardy 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

As part of the National Aerospace Plane Project an analytical model was 
developed to perform calculations for In-line transfer of solid-liquid mixtures 
of hydrogen. This code, called FLUSH, calculates pressure drop and solid frac- 
tion loss for the flow of slush hydrogen through pipe systems. The model 
solves the steady-state, one-dimensional equation of energy to obtain slush 
loss estimates. A description of the code is provided as well as a guide for 
users of the program. Preliminary results are also presented showing the anti- 
cipated degradation of slush hydrogen solid content for various piping systems. 


INTRODUCTION 

Solid-liquid mixtures of hydrogen are currently being considered for use 
in the National Aerospace Plane (NASP) Project. These mixtures, known as slush 
hydrogen, offer the advantages of greater density and heat capacity than normal- 
boiling-point liquid hydrogen. An Important consideration in the use of solid- 
liquid mixtures of hydrogen Is the ability to transfer the fluid through flow 
lines. Previous studies (refs. 1 and 2) have shown that slush hydrogen with a 
solid fraction of 0.5 flows well through lines, valves, and other flow restric- 
tions. However, few studies have shown calculations of the expected solid 
loss due to melting In these flow systems. Ewart and Dergance (ref. 3) per- 
formed such a study for the use of densified cryogenic propellants on advanced 
space transportation systems, but this effort concentrated on line sizes 
greater than 6 in. in diameter and the code it produced is not available in the 
open literature. Therefore, as part of the NASP project, an analytical mode 1 
that performs calculations on the flow of slush in vacuum-jacketed lines was 
developed. This code, called FLUSH (FLow of slUSH), calculates the pressure 
drop and slush hydrogen solid fraction loss for steady-state, one-dimensiona’ 
flow. These calculations can be useful in determining slush transfer parame- 
ters for experimental line analyses as well as parameters important in the 
design of slush transfer lines to the NASP vehicle Itself. 

This report summarizes the method of solution used in the code and pro- 
vides a guide to the users of the code. The program solves the steady-state, 
one-dimensional energy equation and the Bernoulli equation for pipe flow. 

Slush hydrogen and liquid hydrogen properties are obtained by using the GASPl-S 
(ref. 4) physical properties code. The users guide gives input parameters -e- 
quired for FLUSH and example problems on the use of the code. The input re- 
quirements of the program include flow rate, upstream pressure, initial solid 
hydrogen fraction, element heat leak, and element parameters sucn as length 
and diameter. FLUSH also allows the user to specify pipe components such as 
elbows, valves, and straight pipe segments. FLUSH further aPows the user to 
input fr'ction factors and has the capability of calculating heat leak rates 
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for vacuum-jacketed piping. If estimates of the heat leak are not available. 

If standard insulation is desired Instead of vacuum-jacketed piping, the user 
may specify input parameters that will allow the program to calculate heat leak 
for these cases. In addition, the user may vary the flow rate to obtain the 
hydrogen solid fraction loss for a range of flows. A FORTRAN program listing 
of the code is provided, and tables that give basic engineering data for de- 
signing vacuum-jacketed piping systems are presented. 

In addition to the method of solution and the users guide, a summary of 
initial results obtained with the code is presented. The solid fraction degra- 
dation calculated by FLUSH is compared with estimates from previous analyses. 
Tradeoffs between heat leak and friction losses are discussed for various pipe 
sizes and lengths. Finally, a preliminary calculation of the slush loss ex- 
pected at the NASA Lewis Research Center Plum Brook Station's slush hydrogen 
test facility in Sandusky, Ohio, Is presented. 

FLUSH METHOO OF SOLUTION 

The following section describes the method of solution used to calculate 
pressure drop and flow rate for slush hydrogen piping systems. Appendix A de- 
fines the symbols used below. Appendix B gives the actual program listing for 
flush. 


FLUSH solves the steady-state, one-dimensional energy equation given by 
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and the Bernoulli equation given by 
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In the analysis the following assumptions are made: 

1. Steady-state flow exists in the system. 

2. The fluid Is incompressible. 

3. The lines are precooled to triple-point temperature (assumed to be a 
constant) before the slush flows through the lines. 

4. The slush is well mixed before entering the line, and the effects of 
settl ing are neglected. 
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5. The viscosity of slush is the viscosity of liquid hydrogen because the 
viscosity of slush hydrogen is not known. 


6. Heat of fusion it a constant. 

7. The heat leak into the slush is used to melt the solid until all the 
solid is gone; the temperature is assumed to stay constant during the melting 
process. Once all the solid has melted, the heat leak produces a temperature 
rise in the remaining liquid. 

The assumptions that the fluid is well mixed and that the lines are pre- 
cooled are probably the most significant of those listed. 

FLUSH solves the preceding equations for each element specified by the 
user. An element can be a straight pipe length, a valve, a fitting, or any 
other piece of equipment in a piping system. Each element has two nodes; ore 
node is at the beginring of the element, described by a set of initial or 
entrance conditions, and the other node is at the end of the element and is 
characterized by the final or exit conditions. 

The following steps are used in solving these equations for the pressure 
drop through an element: 

1. Initialize temperatures, pressures, slush solid fraction, and heat 
leak, as well as element length, diameter, and height. 

2. Use GASPLUS to obtain initial-element liquid hydrogen properties. 

Pl-Pl<Tl.Pl> <4> 

PI - Ml <Ti ,Pi > (5) 

3. Calculate an effective slush density at the inlet. 

Pmi x * x,P s ♦ (1 ~ *1> p l 5 (6) 

4. Calculate mass flow rate and Reynolds number. The Reynolds number for 
slush is calculated from the liquid viscosities. 
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5. Calculate a friction factor by using either standard liquid friction 
factor correlations (the default selection) or National Institute of Standards 
and Technology (NIST, formerly the National Bureau of Standards) slush correla- 

tions. The standard friction factor correlation used is the Co' brook equation, 
which is a fit of data from the Moody plot of friction factor versus Reynolds 
number . 
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6. Calculate the downstream pressure from Bernoulli's equation. 
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Following calculation of the pressure drop the slush solid fraction exiting the 
element can be calculated. This Is done in the following steps: 

1. Obtain final-element hydrogen properties from GASPLUS. 

PI * Pl<Tf,Pf> (10) 

Ml - Ml (Tf ,Pf ) n]) 

2. Obtain initial- and final-element liquid enthalpies from GASPLUS. 

H 11q,f * H 11q.f(Tf . p f> (12) 

H 1 i q , 1 = H] jg j(Tj ,Pi> (13 > 

3 Calculate initial mixture enthalpy. 

Hmi x. 1 = H 1 i q , i “ x i H f u s v ' 4 

4. Calculate final mixture enthalpy from the energy equation. 


i mi x , f * M mi x , 1 + • 2g c J " g c J 
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The heat leak into the pipe element Q is normally a constant provided by 
manufacturers of cryogenic equipment. The value of W s is zero unless a pump 
is used in the system. 

5. Calculate the change In slush solid fraction due to heat input. 

H 1 iq.f H mix,f , ,& 


6. Calculate the change in sluch solid fraction due to friction losses 


fnct H fus J 

7. Calculate the slush fraction at the element exit. 

xf - *f ,h - Axfrict ( 18 

If slush still exists at the exit, the calculations are eoeated using tne 
solid fraction xf to obtain a new value of p^*. The calculations are 
stopped when convergence Is reached on density and pressure. The converge-:? 
criteria are 0.0001 1 bm/f t 3 and 0.0001 psi, respectively. 

If xf » 0, calculations are performed to find the increase In liquid temce- a 
ture due to heat Input. 
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1. Obtain the Initial element enthalpy, density, and viscosity by using 
GASPLUS. 

2. Obtain the saturated liquid enthalpy and the saturation temperature. 

3. Find the final liquid enthalpy from the energy equation. 

4. Obtain a liquid temperature from GASPLUS. 

T 11q,f - T 11q,f (H liq,f. p f ) <19) 

If T 11q,f > T sat* assume that Tnq t f = T sa t and make a comment in the out 
put. The program does not perform calculations for liquid-gas two-phase flow. 

5. Repeat calculations to obtain convergence on temperature. 

Once the slush solid fraction or the temperature at the output of the element 
has been obtained, the calculations can be repeated for a new element. Follow- 
ing completion of the calculations on a system of elements the initial slush 
flow rate is Increased, If the user chose to vary the flow rate, and the entire 
process is repeated. 

The code provides several options for determining heat transfer into a 
pipe element. If Industry estimates are available, a constant value is used. 

If these estimates are not available for a vacuum-jacketed pipe or if layers of 
insulation are used on the pipe element, the FLUSH code can calculate the heat 
transfer values. 


If the user desires to estimate the heat transfer into a vacuum-jacketed 
pipe element, a subroutine Is provided In the code to perform this calcula- 
tion. The method of calculating this heat transfer rate is discussed by Barron 
(ref. 5) and Is described here. 


The primary source of heat transfer ecross a vacuum-jacketed pipe Is radl 
ation. A cross section of a vacuum-jacket’d pipe Is shown In figure 1. The 
heat rate can be obtained by the equation 


F e F l-2 aA l 


° A l( T 2 - T ?) 


where, for no radiation shields, the emlsslvlty factor F e is defined as 
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Many vacuum-jacketed pipes contain radiation shielding to reduce the radiant 
heat transfer. In this case the emlsslvlty factor Is defined as 


s- . — ♦ — - 1 


(N - 1) 


- 1 ♦ — 

[ c 2 


Ewart and Oergance (ref. 3) point out that actual values of the emlsslvlty fac- 
tor can be considerably higher (a factor of five or greater) than those calcu- 
lated. Therefore. In this study the actual emlsslvlty used will be five times 
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that calculated. For cryogenic vessels the Inner vessel is normally completely 
enclosed by the outer vessel, and Fj _2 =« 1. 

For most vacuum-jacketed lines the vacuum space pressure is lxlO -5 tore or 
less. In this case, as shown by Ludtke and Voth (ref. 6), molecular conduction 
has little effect on the heat transfer. Therefore, calculations for molecular 
conduction In the vacuum space have not been included in the subroutine. 

Note that the Industry estimates of heat leak and those calculated by the 
subroutine are steady-state heat leak values. In an actual transfer line the 
heat leak may be higher owing to chllldown effects. Although this code assumes 
the line Is chilled before the slush Is placed In the line, it Is possible to 
calculate transient effects (ref. 5). This transient analysis, however, is 
beyond the scope of the present work. 

For most ground applications of slush hydrogen, vacuum-jacketed piping 
would be used for transfer. There may be cases, however, where vacuum-jacketed 
line Is not feasible. In this case standard Insulations may be required. Fig- 
ure 2 shows a crosssection of piping with insulation. For conduction through 
insulation the heat transfer Is calculated as 


0 = UAQ(T a j r - Tfi u j(j) 


whe' - e 




Ar a Ar b 

k a A a,lm + k b A b,lm 


Ar, 


k c A c ,1m 


Ar d 

k d A d,lm 


h A 
0 o 


( 24 ) 


where the subscripts a, b, c, and d refer to the pipe material and the 
various types of Insulation. 


The log-mean area Is defined as 



( 25 ) 


c or the case where a pump exists In the flow lines, additional calcula- 
tions are requireo. For a pump the one-dimensional energy equation and 
Bernoulli equation (eqs. (1) and (2)) are solved with the addition of a shaft 
work term as discussed previously. The pump shaft work is defined as 
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60(HP) 

0.02356m 
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In addition, the pressure rise can be calculated once the efficiency of 
the pump and the shaft work are known, as given In equation (2). 


-144AP 

P 


JW s .min 




( 2 /) 


where W s m ^ n Is the minimum shaft work required, or the hydraulic horsepower 
This hydraulic horsepower is defined as 
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w s,min = F/J + W s = nH s 
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FLUSH USERS GUIDE 

The standard Input file for the FLUSH code is divided into two sections 
H 0 n f T,e i 1 f t se f t1 ? n . and t an element description section, in addition, if it is 
desired to calculate the heat leak into any element, a radiation heat transfer- 
parameter section must be added. This section is not required if the heat leak 
rates are known. Also, if standard Insulation or bare pipes are used, two ad- 
ditional sections are required. These are described below. The input file 

"f LUSH -™ p " for the code to run. Additional notes fo? users on 
t e NASA Lewis Research Center VAX system are provided in appendix c. 


In the namelist $PARAMS section the variables are defined as follows: 


FlOMIn 

F LOMAX 
FLOINC 
E 
P 

XI 

ZIN 

NUM 

OUT 

LONG 

XPUT 

FRICT 


minimum flow rate, if it is desired to vary the flow in the anal vis 
gpm v “ 

maximum flow rate, gpm 

flow Increment, gpm 

pipe roughness, ft 

Initial (upstream) pressure, psl 

Initial slush solid fraction 

initial (entrance) height of element 1, ft 

number of elements in system 

output control flag. If OUT < 0, only final output conditions will 
be printed to unit 9; otherwise, output conditions of each eleme-v 
are displayed 


t ->e 


HEAT 

TAMB 

VAC 

XRS 


output control flag for individual variable output. If tONG > o 
Individual variables are printed to unit 8. This is useful in 
examining properties, etc. 

Input control flag. If XPUT > 0, print input data 

NIST friction factor option. If FRICT > 0. use standard friction 
factor correlations; otherwise, use NIST relations for slush 

element heat transfer flag. If HEAT > o. heat transfe- into eacn 
element has units of Btu/hr ft; otherwise, units of Btu/hr a--e used 

ambient temperature, °R 

number of elements that require a calculation for the heat leak rate 
number of radiation shields In vacuum- jacketed line inis 1 < ,-p-, 


for the heat leak calculation. Note 
equation (22) 


XRS is equiva'en 1 


e ; 


to N ’ n 
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NOJAC 


number of elements that do not have vacuum jacketing (see eqs. (23) 
and (24)) 



The element description section of the Input requires the Information nec 
essary to define an element. These variables are defined as follows: 


ELE element number 

TYPE element type. The codes corresponding to the element type are given 
In table 1 

DIA element diameter. In. 

LENGTH element length, ft 

TEMP element temperature, °R 

n element heat transfer rate, Btu/hr or Btu/ft hr 

< flow resistance coefficient, normally used to define valves, pipe 

bends, ana fittings 

HT height of element at Its exit, ft 

Table 2 shows an example of a standard Input file for the FLUSH program^ 

In this example an analysis Is performed on a 100-ft length of 1 .5-in. -diameter 
Schedule 5S vacuum-jacketed pipe. The line Is divided into five sections with 
bellows and bayonet fittings Installed In the line. In addition, two globe 
valves are placed In the line, one near the beginning of the flow and one near 
the end The entire line Is located 25 ft above the ground. It is desired to 
vary the flow from a minimum of 10 to 100 gpm In increments of 5 gpm and find 
how much slush Is lost at the end of the flow length given an initial solid 
fraction of 0.5. In this case conditions at the end of the piping system will 
oe printed for each flow rate; input conditions also will be printed. Standard 
friction factor correlations will be used, and the heat transfer will be de- 
fined in British thermal units per hour foot for each element. The ambient 
temperature is 530 °R, and no heat leak calculations are required ^VAC = O'. 

The variable XRS Is not needed here as no heat leak estimates are being calcu- 
lated, but It Is shown for completeness. 

The element conditions are defined In the second section of the input. 

The temperature of the slush hydrogen Is assumed to be 25. 2 °R. This value 
for the triple-point temperature Is required, rather than the normal value of 
24.84 # R, owing o the limitations of GASPLUS. The heat transfer rates to 
each of the IndU 1 >a1 elements were obtained from estimates made in industiy; 
therefore, the rao it Ion parameter section Is not required. Estimates of neat 
leak to liquid hydrogen for various size pipes are shown in table 3 

As discussed In the literature (ref. 7), valves are normally defined n 
terms of either a flow coefficient C v or a resistance coefficient < where 

c 29.9d 2 1 29 > 

v " V< 
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Typical C v and heat leak values for valves and fittings are given in table 4. 
If < is set equal to zero, the code assumes that a straight vacuum-jacketed 
line is being used. 

Output from this example is shown in table 5. The listing gives the ini- 
tial conditions of each element in the first page of output. The output param- 
eters include the initial volumetric flow rate in gallons per minute, the mass 
flow rate in pounds per second, the final solid fraction, the change In solid 
fraction, and the pressure drop in pounds per square inch. Had all the solid 
hydrogen melted before the piping system exit had been reached, the heat leak- 
age through the lines would be used to increase the temperature of the liquid 
hydrogen remaining. In this case the output under the heading "Change in x" 
would show the final fluid temperature. If the variable OUT had been specified 
as greater than zero, the conditions at each node would have been listed in the 
output. An example of this output is given in table 6. In addition to the 
previously defined parameters this form of output includes the pressure at each 
node (a node is defined as a point at the exit or entrance of an element). 

Table 7 shows an example input file for the case where the heat leak ra*:e 
is calculated for various elements of the previous example. In this case it is 
desired to calculate the rate of heat leak into elements 1 and 9. The variable 
VAC is set equal to 2, and the heat leak rate 0 is set less than zero for 

elements 1 and 9 in the second section of the input. In addition, radiation 

heat transfer parameters must be defined: the element number, the outside 

diameter (o.d.) of the inner pipe, the inside diameter (i.d.) of the outer 
pipe, the emissivlty of the inner pipe ei (El in the input), and the emissivity 
of outer pipe, ti (Eo In the input). The i.d. and o.d. for Schedule 5S stain- 
less steel pipe are given In table 8. For this example a 3-in. pipe is used 
for the outer pipe, and 10 radiation shields between the inner and outer pipes 

are used as defined by XRS = 10. Note that the user must leave a space between 

the element description section and the radiation heat transfer parameter sec- 
tion. Also, the heat leak must be defined in British thermal units per hour 
foot when calculating the heat leak. Typical output for this example is given 
In table 9. 

The output for this example shows a calculated heat leak of 0.552 Btu/hr 
ft, a value somewhat lower than the estimates given In table 3. These differ- 
ences could be the result of differences in the surrounding conditions assumed 
for the Industry estimates. Figure 3 shows the calculated heat leak as a func- 
tion of ambient temperature for various pipe sizes These calculations were 
based on vacuum-jacketed l’ne with 10 radiation shields. Emisslv'ties were as 
sumed to be 0.3 for both the inner and outer lines. This emissivity value cor- 
responds to that of polished stainless steel (ref. 8). Ten radiation shields 
with emisslvltles of 0.03 were assumed for the analysis. Comparing the heat 
leak values given In the figure with Industry estimates of heai leak shows that 
the calculated values are generally between 5 and '3 percent lower than indus- 
try estimates, assuming an ambient temperature of 340 J R. Therefore, the cal- 
culated values of heat leak are expected to be dose to those given by industry 
as design values. 

When standard Insulation is used, the type of insulation and the thickness 
must be specified. An example of such a case is shown in table 10. Here, in 
the first part of the Insulation parameters, ElE is the element number, o.D. is 
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the pipe outer diameter, and T1 , T2, and T3 are the thicknesses of the insula- 
tion. Up to three Insulation layers are allowed by the code. In the second 
section of the Insulation parameters the material types are specified. The 
material types corresponding to the codes are given In table 11. In this 
example five elements have either standard Insulation or bare pipes, as speci- 
fied by N0JAC»5. Various combinations of pipe materials and insulation schemes 
are shown. Note that the value of 0 given In the element description section 
will be replaced by that calculated In this section; therefore, the 0 speci- 
fied will not be that given in the output. The output for this example is 
shown in table 12. 

Finally, an example is presented in table 13 for the case where a pump 
exists in the flow system. In order to specify a pump in the line, the element 
type is set as either PP or PH. Element type PP signifies a pump with pressure 
rise and efficiency specified, while PH represents a pump with brake horsepower 
and efficiency given. In the example shown element 2 is set to type PH. In 
this case the horsepower is defined in the column under Q and the efficiency 
is input in the column under K. Also, element 4 Is specified as type PP. The 
pressure drop for the pump in element 4 is input under the 0 column and effi- 
ciency is placed in the column under K. Table 14 shows the output for the 
example; slush loss would occur in the pump owing to energy input by the pump. 

The user should be aware that there are some limits on the code. These 
’imits are summarized as follows: 

1. The triple-point temperature of slush must be set equal to 25.2 °R, 

"at h er than 24.8 °R, owing to limitations on GASPLUS. 

2. The exit pressure on the piping system must be greater than 1.1 psi 
(approximately the triple-point pressure of hydrogen). If this limit is not 
met because the inlet pressure was chosen to be too low, a message appears in 
the output file. 

3. The maximum number of pipe elements the user can input Is 90. 

4. When all the slush has melted and the heat leak into the system causes 
a rise In temperature, the code Is limited to a final temperature equal to that 
of the normal-boiling-point saturated liquid; gas-liquid two-phase flow is not 
included In FLUSH. 

5. When using the NIST friction factor correlations, a maximum Reynolds 
number of only 1x10° Is possible because of limits on the data. 

6. Up to three layers of insulation can be used for the pipe if vacuum 
jacketed piping Is not used. 


RESULTS 

In order to obtain confidence in the code, a comparison was made w : tn >e 
suits from a previous analysis performed by Ewart and Dergance (ref. 3). In 
this study propellant heating due to heat leak and friction losses was consid- 
ered, as In the present study. The losses due to friction were assumed to be 
the same as those for liquid hydrogen In the Ewart and Dergance study; there- 
fore, standard friction factor correlations were used in this comparison. The 
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line considered was a 12-1n. -diameter Schedule 5S pipe, 1750 ft long, with a 
constant heat leak of 12 Btu/hr ft. This heat leak value was calculated for a 
vacuum-jacketed pipe with multilayer Insulation. 

Figure 4 compares the slush loss calculated by FLUSH with that of the 
Ewart and Dergance study. As can be seen from the data, the results are nearly 
the same over the entire range of flow rates. Also, as discussed In the previ- 
ous study, a minimum slush loss is reached at a flow rate between 3000 and 
4000 gpm owing tradeoffs in environmental heating and friction heating. At 
lower flow rates heat leak Into the line is the major source of slush loss. 

At higher flow rates friction heating becomes significant. For example, in the 
Ewart and Dergance comparison at 500 gpm, the solid fraction loss due to envi- 
ronmental heating Is 0.043, while that due to f r1ction losses is essentially 
zero. At 10 000 gpm, the solid fraction loss due to environmental heating is 
0.013, while the loss due to friction is 0.012. For cases with less heat leak 
the slush loss due to friction becomes the dominant loss mechanism. 

Following a comparison against previous results the FLUSH code was run 
using various line sizes and pipe lengths. The results of these runs are shown 
in figure 5 for line lengths of 100, 200 and 500 ft, respectively. Each line 
was assumed to be only a straight section of Schedule 5S pipe, with no valves 
or fittings. As seen in each of these figures, a minimum solids loss was 
reached as in the case discussed previously. As the line size is decreased, it 
appears that the minimum occurs at a lower flow rate and at a higher solid 
fraction loss. This Is due to the Increased friction loss for smaller pipes, as 
can be seen in the equation 


whei e F is the energy loss due to friction. The effect of a smaller pipe 
cross-sectional area is to increase the velocity which, coupled with a de- 
creased diameter, leads to an increased friction loss. Hence, for a small pipe 
size the friction loss will begin manifesting Itself at a lower flow rate, and, 
as the flow rate increases, the smaller pipe size will present higher solid 
fraction losses than do the larger pipe sizes. At low flow rates the larger 
pipes exhibit large slush losses. This is due to the higher surface area and 
hence higher heat leak Into the larger pipes. For longer pipes the solid frac- 
tion losses Increase owing to Increases in both the friction loss and the envi- 
ronmental loss. 

A piping system Is presently being designed for the transfer of slush 
hydrogen at NASA Lewis Research Center's Plum Brook Station (ref. 9). This 
system, located in the K-SIte facility, will require the transfer of slush hy- 
drogen from a slush generator located approximately 100 ft r rom the slush test 
tank through 1 .5-1n. -diameter Schedule 5S pipe. An analysis of this piping 
system was performed by using the FLUSH code. The current system design con- 
tains two valves and five 20-ft sections of pipe. Bayonet f ‘ittings are used 
to connect these lines, and bellows are deluded in each of the lines. The 
heat leak used for the piping was obtained from industry estimates, provided 
in tables 3 and 4. The heat leak to the Pel lows was assumed to ce tnat of t l 'e 
flexline. The roughness of the pipe was assumed to be 0.00015, typically the 
roughness of commercial steel. The flow rate was varied between 10 and 1 30 gpm 
and values of the slush loss were obtained. 







The results of the slush loss analysis for K-Slte are shown in figure 6. 
The results show that a fairly broad minimum occurs for this system. Between 
35 and 60 gpm the solid loss is approximately 0.5 percentage point. In other 
words, the expected slush quality, if the fluid is transferred it a rate 
between 35 and 60 gpm, will be 49.5 percent when it reaches the test tank if 
the slush 1s Initially 50 percent solid hydrogen. This represents a small 
loss when using standard vacuum-jacketed pipe with multilayer insulation. 
Therefore, helium transfer technology is not required for this system. 

Finally, a comparison was done between the K-Slte system for standard 
friction factor correlations and the same system for NIST friction factor re- 
sults (refs. 1,2, and 10). The experimental NIST friction factor results, 
shown in figure 7, were obtained on a transfer line 0.652 in. i.d., and 
48.4 ft long. The data were obtained for Reynolds numbers less than 8x10^ and 
extrapolated to the lower flow rates. The solid line in the figure represents 
a standard friction factor correlation used by NIST. Because these friction 
factors were obtained empirically, it is not known how the friction factors 
will change outside the range of Reynolds numbers tested. In addition, the 
NIST results were obtained for smooth pipes. Therefore, the <-SUe case dis- 
cussed above was run again assuming a smooth pipe. The comparison is shown >n 
figure 8 for an initial slush solid fraction of 0.5. As expected, the slush 
loss calculations are nearly identical at flow rates less than the minimum 
value because environmental heating dominates the loss of slush. At higher 
flow rates the slush loss estimates are quite close owing to the s i mi la- 
values of the friction factor. For the case using the NIST correlations 
values of slush loss are not shown for flow rates greater than '75 gpm owing 
to the constraint that the Reynolds number must be less than 1x10®. 

An additional issue that has not been addressed is that of critical veloc- 
ity, the velocity where the solid hydrogen particles begin to settle, a previ- 
ous study (re f . 1) has shown that for a 0. 652-in. -i .d. pipe the critical 
velocity is 1.5 ft/sec. If this value is applied to larger pi^e sizes, the 
c i t i ca 1 flow rates range from 5.2 gpm for a 1-in. -diameter Sci.edu 1 e 5S pipe to 
568 gpm for a 12-1n. -diameter pipe, as shown in table 15. Below the critical 
flow rate the fluid becomes stratified, and the assumptions of one-d'mensional 
flow and a well-mixed fluid are no longer valid. Therefore, the values pro- 
duced by the code are not expected to be accurate at flow rates below the crit- 
ical value. It is clear that more testing is required to obta'n the critical 
velocities at larger pipe sizes as the values in table 15 are extrapolated from 
data for small lines. It has been postulated that the critical velocity de- 
pends on pipe diameter with the dependence ranging from direct proportionality 
to diameter to the 0.15 power (ref. 11). 


CONCLUDING REMARKS 

An analytical model was developed as part of the National Aerospace Plane 
Project to predict pressure drop and hydrogen solid fraction loss for a range 
of pipe sizes and flow rates. This code, called FLUSH, has di-ect aDDiicat’cn , 

to the design of piping systems for slusn hydrogen exper imenta ’ test systems | 

and for the flight and ground handling systems of the NASP vehicle. r he code 
solves the one-dimensional , steady-state heat tr ansfer eauation *o f ; n.j the 
solid fraction loss. A listing of the program is provided, anj a guide to tne 
code's use Is presented. 
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Results using FLUSH show that this code agrees well with previous analy- 
ses. Also, the data show that at low flow rates heat leak effects dominate the 
solid fraction loss, while at higher flow rates friction heating causes most 
of the solid hydrogen degradation. These two effects result In a minimum slush 
solid fraction loss for each pipe size. Further, as the pipe diameter Increases, 
this minimum slush solid fraction loss decreases owing to the decreasing fric- 
tion losses. An analysis of the NASA Plum Brook slush facility transfer line 
showed small losses for the piping system design being considered. The analy- 
sis will be compared with the test results upon completion of the facility. 
Howev?r, solid hydrogen losses appear to be In an acceptable range for fully 
cnllled vacuum-jacketed lines designed for liquid hydrogen technology. 

Finally, comparing slush loss with standard friction factor correlations and 
with friction factor relations obtained In previous experiments showed similar 
values over the range of flow rates in this study. 
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APPENDIX A 



SYMBOLS 


A a , 1 m \ 



A b,lm | 

log mean area of material a, b, c, or d, ft 2 


A c,lm ^ 



Ad , 1m 



A i 

inside surface area of pipe, ft 2 


A o 

outside surface area of pipe, ft 2 


A 1 

area of surface 1 per foot, ft 2 /ft 


a 2 

area of surface 2 per foot, ft 2 /ft 


Cv 

flow coefficient 


D 

pipe diameter, ft 


d 

pipe diameter. In. 


F 

energy loss due to friction, ft lbf/lbm 


Fa 

emissivity factor 


F i-2 

configuration factor 


f 

friction factor 


g 

acceleration due to gravity, 32.2 ft/sec 2 


9c 

32.2 lbm ft/lbf sec 2 


AH 

change in enthalpy, Btu/lb 


HP 

brake horsepower, hp 


HfuS 

latent heat of fusion, Btu/lb 


H 11q,f 

enthalpy of liquid exiting element, Btu/lb 


Hljq.l 

enthalpy of liquid entering element, Btu/lb 


H m1 x ,f 

enthalpy of solid-liquid mixture exiting element. Btu/lo 


H m1 x , 1 

enthalpy of solid-liquid mixture exiting element. Btu/lb 


h i 

Inside heat transfer coefficient, Btu/hr ft 2 °R 


n 0 

outside heat transfer coefficient, Btu/hr ft 2 9 R 
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! H 


l 4 


J 


mechanical equivalent of heat, 778.2 ft-lbf/Btu 
resistance coefficient 



m 


N 


A? 

Pf 

Pi 

0 

q 

Re 
A’ a 

A'o 

Ar c 

Ar 


thermal conductivity of material a. 


b, c, or d, Btu/hr ft R 


element length, ft 
mass flow rate, lb/hr 
number of radiation shields 
change In pressure, psi 
pressure at element exit, psi 
pressure at element entrance, psi 
heat leak Into element, Btu/hr 
volumetric flow rate, gpm 
Reynolds number (also R> 
thickness of material a (r a - rp, ft 

thickness of material b (rjj - r a ), ft 

thickness of material c (r c - rj>), ft 

thickness of material d (rj - r c ), ft 


r a* r b 1 

outside radius of material a, b, c, or d, ft 


r c* r d ) 



n 

Inside radius for log-mean radius calculation, 

ft 

r 2 

outside raldus for log-mean radius calcuatlon, 

ft 

^air 

air temperaure, °R 


T f 

temperature at element exit, ®R 


i u * a 

fluid temperature, ®R 


T i 

temperature at element entrance, °R 











1 

1 


i 
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T 11q.f 

Tsat 

Tl 

^2 


U 

V 

w s 

w s ,mln 
Xf 

V ,h 


* i 


.}Z 

c s 

ci 

c 2 

n 

M 

Ml 

P 

Pmi « 


PI 

<7 


liquid temperature at element exit, °R 
saturated liquid temperature, °R 

temperature of surface 1 (equal to fluid temperature', a R 
ambient temperature, *R 

overall heat transfer coefficient, Btu/hr ft 2 °R 
fluid velocity, ft/sec 
shaft work, Btu/lb 

minimum shaft work (hydraulic horsepower), Btu/lb 
solid fraction at element exit 

solid fraction at exit (heat leak considerations only' 

change in solid fraction due to friction losses 

solid fraction at element Inlet 

change In element height, ft 

emlsslvlty of radiation shielding 

emlsslvlty of surface 1 

emlsslvlty of surface 2 

pump efficiency 

triple-point hydrogen viscosity 
liquid viscosity, ibm/ft sec 
slush density, lbm/ft^ 
solid-liquid mixture density, ibm/ft^ 
sol Id density, ibm/ft^ 
liquid density, ibm/ft^ 

Stefan-boltzmann constant, 0.173x10*® Btu/hr ft' -'P J 




16 





REFERENCES 

1. Slndt, C.F.; and Ludtke, P.R.: Slush Hydrogen Flow Characteristics and \ 

Solid Fraction Upgrading. Advances In Cryogenic Engineering, Vol . 15., ; 

K.D. Tlmmerhaus, ed.. Plenum Press, New York, 1970, pp. 382-390. • 

2 . Slndt, C.F.; Ludtke, P.R.; and Daney, O.E.: Slush Hydrogen Fluid Charac- i 

terlzatlon and Instrumentation. NBS Technical Note 377, Feb. 1969. 3 

■j 

3. Ewart. R.O.; and Dergance, R.H.: Cryogenic Propellant Denslf icatlon j 

CR U ?59438 MC ?978~ 586 ’ Maft1n Mar,etta Cr °Pi NASA Contract NAS3-2101 4) NASA \ 

i 

4. Fowler, J.R. : GASPLUS User's Manual NASP CR- 1012, Mar. 1988. ] 

5. Barron, R.F.: Cryogenic Systems. 2nd ed., Oxford University Press Ox- i 

ford, 1985. 


5. Ludtke, P.R.; and Voth, R.O.: A Study of LC-39 Cryogenic Systems. (NBS 

Report 10705. National Bureau of Standards, NASA CR- 1 26386 , 1971. 

1 . Flow of Fluids Through Valves, Fittings, and Pipe. Crane Co., Hydro-A’-e 
Division, Chicago, IL, TP-410, 1965. 

3. Siegel, R.; and Howell, J.R.: Thermal Radiation Heat Transfer. 2nd ed.. 

McGraw-Hill, New York, 1981. 

9. Hannum, N.P.: Technology Issues Associated with Fueling The National Aero- 

space Plane With Slush Hydrogen. NASA TM- 101 386 , 1988. 

19. Slndt, C.: A Summary of the Characterization Study of Slush Hydrogen. 

Cryogenics, vo’ . 1C, no. 5. Oct., 1970, pp. 372-380. 

! ‘ ■ Voth - R -0-. et al: On Process Design Criteria, Data and Methodology Relat- 

ed to Slush Hydrogen Technology. NASP TM-1006, 1987. 


I 

* 

J 



17 





APPENDIX B 
PROGRAM LISTING 


18 





o o o o 



PROGRAM FLUSH 

C 

c 

C THIS PROGRAM CALCULATES PRESSURE LOSS AND SLUSH HYDROGEN 

C SOLID FRACTION LOSS IN STEADY STATE FLOWS IN PIPE SYSTEMS 

C (SIMPLE PUMP EQUATIONS ARE INCLUDED) 

C 

C THE PROGRAM SOLVES THE STEADY-STATE, ONE DIMENSIONAL 

C ENERGY EQUATION. 

C 

C H2-H1 * (V2**2-V1**2) /2GC ♦ (G./GC) (Z2-Z1) ■= Q - WS 

C 

(’ AND THE BERNOULLI EQUATION 

C 

C (P2 Pl)/RHO * (V2**2 Vl*»2)/2GC ♦ (G/GC) (Z2-Z1) » F ♦ WS - 0 

C 

C 

C 

C 

C TERRY HARDY 

C NASA LEWIS RESEARCH CENTER 

C CLEVELAND, OHIO 44135 

C (216) 433-2411 

C 

C 4-17-89 

C 

C- 

c 

DATA ID/ ’PH2’/ 

DIMENSION EDIA(90) ,ELBN(90) ,ET(90) ,EX1 (90) ,BZ2(90) ,EQIN(90) , 

• ERHO(QO) ,EP1 (90) ,EMD0T(9O) ,BQUAL(90) ,EDP2(90) , 

* N0DB(9O) ,TS(90) ,HF1 (90) ,ETNBW(90) ,ECLD(90) ,RE1(90) , 

♦ NVAC(90) ,BVDI(90) ,BVD0(9O) ,BEI(90) ,EE0(9O) ,PTYPA(90) , 

• QFL(90) ,QL(90) ,PXTYP(10,90) ,INSUL(90) 

COMMON /FLUX/ EP0D(9O) ,ETI1 (90) ,ETI2(90) ,BTI3(90) , 

NNOVJ (90) , NO JAC , NXMAT (20 , 90) , 

. EXLEN (90) , EXDIA (90) , EXT (90) , TAMB , P , X 1 , FLOV , EPID (90) , K 

C 

REAL MDOT , LEN , MPHR , LONG 
CHARACTER* 12 PXTYP 
CHARACTER* 10 TIM, TONE, TAM, TANE 
CHARACTER*4 PTYP , TYPE , PTYPA 


READ IN INPUT DATA 


C 

OPEN (UNIT=7 , FILE- ’ FLUSH. INP STATUS= ’ OLD ’ ) 

C 

C PARAMS NAMELIST SECTION 

C 

NAMELIST / PARAMS / PL0MIN,FL0MAX,FL0INC,E,P,X1 .ZIN.NUM.OUT, 
* LONG, XPUT.FRICT, HEAT, TAMB, VAC, XRS, NO JAC 
READ (7, PARAMS) 

READ (7, 800) 

POO FORMAT (A80) 

NUMB = NUM 
TOTLEN - 0. 
c 
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■**»> «. wt*: * '* ' * - . 


■f 


c 

c 


r 


I 


300 


I 


I 

S- 


c 

c 

c 

c 


31 


C 

r 

(' 

c 


c 

(' 

c 

c 

c 

c 

c 

c 

0 

c 

c 

c 

c 

c 

c 

c 

u 


r 

r 


ELEMENT DESCRIPTION SECTION 
DO 31 1=1, NUMB 

READ (7 ,300) NELEM , TYPE , DIAM , LEN , T , QIN , CLD , ZHT 
FORMAT (14 , A4 , 6P10 . 2) 

NODE (I) = NELEM 
PTYPA(I)= TYPE 
EDIA(I) = DIAM 
ECLD(I) = CLD 
ELEN (I) = LEN 
ET(I) = T 
ETNEW(I)= T 
EQIN(I) = QIN 
EZ2(I+1)= ZHT 
EZ2(1) = ZIN 

EX1(I) = XI 
TOTLEN = LEN ♦ TOTLEN 
CONTINUE 


FLOMIN= MINIMUM FLOW RATE, GAL/MIN 

PLOMAX= MAXIMUM FLOW RATE, GAL/MIN 

PLOINC= INCREMENT POR FLOW INCREASE, GAL/MIN 
QIN = CONSTANT HEAT LEAK TO VJ PIPE,BTU/HR OR HTlUlii FT 
IF A PUMP ELEMENT, PH OR PP, IS USED 
QIN = BRAKE HORSEPOWER, HP (PH) OR HEAD LOSS, PSI (H’t 
TYPE = ELEMENT TYPE 

PI PIPE ELEMENT 

PP PUMP, PRESSURE LOSS SPECIFIED 

PH PUMP, BRAKE HORSEPOWER SPECIFIED 

BE BELLOWS 

VA STANDARD GLOBE VALVE 

BA BAYONET FITTING 

EL DO DEGREE ELBOW 

B4 45 DEGREE ELBOW 

FX FLEXIBLE TUBING 

YV Y -PATTERN GLOBE VALVE 

AV ANGLE VALVE 

GV GATE VALVE 

BV BALL VALVE 

E = PIPB ROUGHNESS 

CLD = RESISTANCE COEFFICIENT, K (VALVES, ELBOWS , ETC ) 

K = F (L/D) 

CLD = 0 => STRAIGHT PIPE 

IF PUMP BLRMBNT (PP OR PH), CLD = PUMP EFFICIENCY 
DIAM = PIPB DIAMBTBR, INCHES 

T = INITIAL FLUID TEMPERATURE, R 

P = INITIAL (UPSTREAM) PRESSURE, PS1A 

XI = INITIAL SLUSH SOLID FRACTION 

LEN = PIPB LENGTH, FT 

NUM = NUMBER OF ELEMENTS 

TOTLEN = TOTAL PIPE LBNGTH, FT 
ZHT = FLUID HEIGHT AT NODE N, FT 

ZIN = FLUID HEIGHT AT NODE 1, FT 

OUT = OUTPUT CONTROL FLAG (UNIT 0) 

< 0 => SHORT OUTPUT (FINAL CONDITIONS ONLY i 

LONG OUTPUT CONTROL FLAG FOR INDIVIDUAL VARIABLE. OUTPUT 
(UNIT 8) 

< 0 => NO INDIVIDUAL VARIABLE OUTPUT 
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c 

c 

c 

c 

c 

c 

0 

c 

c 

c 

G 

c 

c 

C-- 

c 

c 

c 

c 


801 

802 


14 

301 


13 

95 


XPUT = INPUT PRINT CONTROL FLAG 

< 0 => DO NOT PRINT INPUT DATA 
HEAT = HEAT TRANSFER PER ELEMENT FLAG 

< 0 => BTU/HR, > 0 => BTU/HR-FT 
FRICT = NBS SLUSH FRICTION FACTOR OPTION 

< 0 => USE NBS CORRELATION FOR F 

VAC = NUMBER OF ELEMENTS WHICH REQUIRE A CALCULATION FOR 
THE HEAT LEAK RATE (VACUUM JACKETED PIPING) 

XRS = NUMBER OF RADIATION SHIELDS FOR VJ LINE CALCULATIONS 
TAMB = AMBIENT TEMPERATURE, R 

NOJAC = NUMBER OF ELEMENTS WHICH REQUIRE A CALCULATION FOR 

THE HEAT LEAK RATE (BARE PIPE OR STANDARD INSULATION) 


CALCULATE THE HEAT LEAK RATE THROUGH VACUUM- JACKETED PIPES 
IF VAC > 0 

jg _ j 

IF (VAC .LE. 0.0) GO TO 95 
DO 802 1=1,4 
READ(7,801) 

FORMAT (A80) 

CONTINUE 

DO 13 IB=1,NUMB 

NEWQ = EQIN(IB) 

TINRAD = ET(IB) 

IF (NEWQ .LT. 0.0) GO TO 14 
GO TO 13 

READ (7, 301) NELEM, DIN, DOUT, El, EO 
FORMAT(I4,2X,4F10. 3) 

NVAC(IE) = NELEM 
EVDI(IE) = DIN 
BVDO(IE) = DOUT 
EEI(IE) = El 
EEO(IE) = EO 

CALL VACUUM (QIN 1 , TAMB , XRS , D IN , DOUT , TINRAD , El , EO ) 

EQIN(IB) = QINl 
IE = IE + 1 
CONTINUE 
CONTINUE 


CALCULATE HEAT LEAK FOR PIPE WITHOUT VACUUM JACKET IF 
NOJAC > 0 


IF (NOJAC .LE. 0) GO TO 160 

DO 311 1=1,4 

READ(7,310) 

310 FORMAT (A80) 

311 CONTINUE 

DO 213 1=1, NOJAC 

READ (7,312) NELB, POD, TI1,TI2,TI3 

312 FORMAT(I4,4F10.2) 

N = NELE 

EPOD(N) = POD/ 12.0 
ETI1(N) = Til 
ETI2(N) = TI2 
ETI3(N) = TI3 
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NNOVJ(I)= NELE 

213 CONTINUE 

DO 214 1=1,3 
READ(7,314) 

314 FORMAT (A80) 

214 CONTINUE 

DO 215 I=l,NOJAC 

READ (7 ,315) NELE, MAT1, MAT2, UAT3, MAT4 

315 FORMAT (14 , 6X , 15 , 5X , 15 , 5X , 15 , 5X , 15 , 5X) 

N = NELE 

NXMAT(l.N) = MAT1 
NXMAT(2,N) = MAT2 
NXMAT(3,N) = MAT3 
NXMAT(4,N) = MAT4 
NNOVJ(I) = NELE 
EXLEN(N) = ELEN(N) 

EXDIA(N) = EDIA(N)/12.0 
EPID(N) = EDIA(N) /12.0 
EXT(N) = ET(N) 

215 CONTINUE 
FLOV = FLOMIN 

DO 216 I=l,NOJAC 
K = NNOVJ(I) 

CALL QFLUX(QFL(I) ,QL(I)) 

IF (HEAT .LT. 0.) GO TO 217 
GO TO 218 

217 EQIN(NNOVJ(I)) = QFL(I) 

GO TO 210 

218 EQIN(NNOVJ(I)) = QL(I) 

219 CONTINUE 

216 CONTINUE 

150 CONTINUE 

C 

C 

WRITE (0,500) 

500 FORMAT (’ FLUSH CODE’/) 

C 

C 

C CALL SYSTEM DATE AND TIME ROUTINES 

C - 

C 

CALL DATE (TIM) 

CALL TIME (TONE) 

WRITE (9, 130) TIM, TONE 
130 FORMAT (/’ ’, A10,2X,A10) 


C WRITE INPUT PARAMETERS TO OUTPUT FILES IF XPUT 0 

C, .. . 

c 

IF (XPUT .LT. 0.0) GO TO 42 
WRITE (9, 406) 

405 FORMAT (/’***************•**************»**♦♦*** * ♦ *♦ • • • 

***** * ♦*•**....**’/) 

WRITB(9,400) 

400 FORMAT (/’ INPUT PARAMETERS ’) 

WRITB (0,401) 

401 PORMAT(/’ ELEMENT’, 2X, ’TYPE’, 2X,’D1A. , IN. ’ ,3X, ’LENGTH. FT 
* ,3X, ’TEMP. , R’, 4X, ’ Q, BTU/HR-FT’ ,3X, ’ K F L/D ’/) 

DO 43 IA=1 ,NUM 
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433 

434 

435 


402 

403 

43 

C 

c 

c 

410 

406 


407 
06 

408 


C 

c 

(1 


97 

411 

412 

413 


414 
106 

419 

415 


IF (HEAT .LT. 0.0) GO TO 434 
GO TO 433 
QALT = EQIN(IA) 

GO TO 435 

QALT = EQIN (IA) /ELEN (I A) 

CONTINUE 

WRITE (9, 402) NODE(IA) ,PTYPA(IA) ,EDIA(IA) , ELEN (I A) ,ET(IA) , 
QALT ECLD(IA) 

F0RMAT(I4,5X, A4.3X.F7. 3,5X,F8.3,5X ,F8.2 .4X.F10. 3.4X.F10. 3) 
WRITE (9, 403) 

FORMAT ( ’ 

CONTINUE 

IF VAC > 0 PRINT RADIATION HEAT TRANSFER PARAMETERS 
IF (VAC .LE. 0) GO TO 97 

WRITE (9,410) mmni>n „ 

FORMAT (/’ RADIATION HEAT TRANSFER PARAMETERS ’) 
um jTR (Q 4061 

FO RMA T (/ * ELEMENT’ ,4X, ’INNER O.D. ’ ,2X, ’OUTER I.D.’,6X, 

’ El \8X,’E0’) 

DO 96 j C-~ 1 VAC 

WRITE (9, 407) NVAC(IC) ,EVDI(IC) .EVDO(IC) ,EEI(IC) ,EEO(IC) 
FORMAT (/I4 , 7X , F8 . 3 , 5X , F8 . 3 , 5X , F8 . 5 , 3X , F8 . 5) 

CONTINUE 
WRITE (9, 408) 

FORMAT ( ’ - 


■’) 


971 


IF NOJAC > 0 PRINT BARE PIPE/STANDARD INSULATION PARAMETERS 
IF (NOJAC .LE. 0) GO TO 197 

FORMAT (/^' BARE PIPE OR STANDARD INSULATION PARAMETERS’) 

FORMAT (/^ELEMENT ’ ,4X, ’ PIPE 0 D . ’, 3X INSULATION 1’, 

4X, -INSULATION 2’ ,4X, ’INSULATION 3’) 

WRITE (0 413) 

FORMAT (26X , ’THICKNESS , FT ’ , 3X , ’ THICKNESS , FT ’ , 3X , 
’THICKNESS, FT’) 

DO 190, 1=1, NOJAC 
K = NNOVJ(I) 

FORD = EP0D(K)*12. 

WRITE (9, 414) K, FORD.ETIl(K) ,BTI2(K) ,ETI3(K) 
FORMAT(/I4,10X,F7.4,6X,F7.4,9X,F7.4,9X,F7.4) 

CONTINUE 

WRITE(9,419) 

FORMAT (A80) 

FOR^T( /^ELEMENT’ ,4X, ’PIPE MATERIAL’ ,4X, ’INSULATION 1*. 
4X, ’INSULATION 2’, 4X , ’ INSULATION 3’) 

DO 954 1 = 1 .NOJAC 
K = NNOVJ(I) 

IF (NXMAT(2,K) .LE. 0.0) GO TO 971 
IF (NXMAT(3,K) .LE. 0.0) GO TO 972 
IF (NXMAT(4,K) .LE. 0.0) GO TO 973 
GO TO 974 
INSUL(I) = 1 
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GO TO 975 

972 INSUL(I) = 2 

GO TO 975 

973 INSUL(I) = 3 

GO TO 975 

974 INSUL(I) = 4 

975 CONTINUE 

954 CONTINUE 

C 

DO 950 I = l.NOJAC 
DO 951 J = 1,4 
K - NNOVJ(I) 

r 

C MATCH MATERIAL CODE TO WORD DESCRIPTION 

C 


IF (NXMAT(J,K) .EQ. 101) GO TO 911 
IF (NXMAT(J.K) .EQ. 102) GO TO 912 
IF (NXMAT(J.K) .EQ. 103) GO TO 913 
IF (NXMAT(J.K) .EQ. 104) GO TO 914 
IF (NXMAT(J.K) .EQ. 201) GO TO 921 
IF (NXMAT(J.K) .EQ. 202) GO TO 922 
IF (NXMAT(J.K) .EQ. 203) GO TO 923 
IF (NXMAT(J,K) .EQ. 204) GO TO 924 
IF (NXMAT(J.K) .EQ. 206) GO TO 925 
IF (NXMAT(J.K) .EQ. 301) GO TO 931 
IF (NXMAT(J.K) .BQ. 302) GO TO 932 




GO 

TO 942 

911 

PXTYP (J,K) 
GO TO 970 

— 

’TEFLON’ 

912 

PXTYP (J,K) 
CO TO 970 

= 

’304 ST STEEL’ 

9 1 3 

PXTYP (J.K) 
GO TO 970 

= 

’CARBON STEEL’ 

914 

PXTYP (J.K) 
GO TO 970 


’6063 ALUM. ’ 

9 2 1 

PXTYP(J,K) 
GO TO 970 


’FOAM GLASS’ 

9 2 2 

PXTYP (J.K) 
GO TO 970 

— 

’SILICA’ 

923 

PXTYP (J.K) 
GO TO 970 

= 

’POLYURETHANE’ 

924 

PXTYP (J,K) 
GO TO 970 

= 

’POLYSTYRENE’ 

92b 

PXTYP (J.K) 
GO TO 970 

= 

’FOAM RUBBER’ 

931 

PXTYP (J.K) 
GO TO 970 

= 

’FIBER GLASS’ 

932 

PXTYP (J.K) 
GO TO 970 

= 

’PERLITE’ 

942 

PXTYP (J,K) 

- 

1 > 

970 

CONTINUE 



951 

CONTINUE 



950 

CONTINUE 




DO 426 1 = 1, NO.) AC 
K NNOVJ(I) 

WRITE ( 9 , 427) K , PXTYP ( 1 , K) , PXTYP (2 , K ) , PXTYP (3 ,K) , PXTYP(4 . K ) 
FORMAT (/I4,7X,A13,4X,A13,3X,A13,3X,A13) 

476 CONTINUE 
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WRITE (9, 428) 
428 FORMAT (’ 


197 WRITE (9, 404) 

404 FORMAT (/ ’********»****♦******»»*****,,**,**,********,,, 
42 CONTINUE /,) 


FLOW - FLOMIN 
TO - T 
PIN1 - P 
XIN = XI 

C 

IF (OUT .GE. 0.0) GO TO 65 
WRITE (9, 522) 

522 FORMAT (’1’,’ FLUSH CODE OUTPUT’) 

WRITE (9, 527) 

527 FORMAT (/’ 


WRITE(9,524) DIAM.FLOW 

524 FORMAT(/’ PIPE DIAMETER = ’ ,F6.2, IX, ’IN. ’ ,9X, ’FLOW 
♦ F7 . 2 , 2X , ’ GPM ’ ) 

WRITE(9,52S) TOTLEN.XIN 

525 FORMAT(/’ TOTAL PIPE LENGTH - ’ .F8.2.1X, ’FT’ 4X ’INITIAL 

* SOLID FRACTION = \F5.3) INI UAL 

WRITE (9, 529) E, TO 

>29 FORMAT(/’ PIPE ROUGHNESS - ’ ,F9. 6, 9X, ’INITIAL TEMPERATURE 

* F7 . 2, IX, ’ R ’ ) 

WRITE (9 ,537) 

537 FORMAT (/’ 


IF (XIN .LE. 0.0) GO TO 66 
WRITE(9,523) 

£23 , FL0W ’ gp “’.3X,’FL0W, LB/S’, 6X,’ FINAL SOLID’, 5X. 

♦ CHANGE IN X ’ ,6X, ’PRBSSURE ’) 

WRITE (9, 528) 

528 P0RMAT(32X, ’FRACTION’, 27X ’DROP, PSI’/) 

GO TO 65 

C 

66 WRITE (9, 533) 

533 , PL0W ’ gp “’.3X,’FL0W, LB/S’, 6X,’ FINAL SOLID’, 5X. 

* PINAL TEMP, R ’ ,4X, ’PRESSURE ’) 

WRITB(9,538) 

538 F0RMAT(32X, ’FRACTION’ ,27X, ’DROP, PSI ’ /) 

65 CONTINUE 


50 CONTINUE 

C 

C RERUN CALCULATIONS FOR A NEW FLOW RATE 

C 

r 

c 

C CALCULATE NEW HEAT LOSS DUE TO CONDUCTION FOR PIPES 

C WITHOUT VACUUM JACKETING 
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C — 

c 

IF (NOJAC .LE. 0) GO TO 151 
DO 220 1=1, NOJAC 
FLOV = FLOW 
K - NNOVJ(I‘t 

CALL QFLUX(QFL(I) , QL (I ) ) 

IF (HEAT .LT. 0.) GO TO 221 
GO TO 222 

221 EQIN(NNOVJ(I) ) = QFL(I) 

GO TO 223 

222 EQ1N(NN0VJ(I)) * QL ( I ) 

223 CONTINUE 

C WRITE(15,271) FLOW, K, EQIN(K) , FLOV 

C 271 FORMAT(F10. 3 ,3X, 14 ,5X ,F14 . 3,5X,K14 . 3) 

220 CONTINUE 

151 CONTINUE 

C 

c 

C RESET TEMPERATURES TO INITIAL VALUES FOR NEW FLOW RATE 

C 

0 

DO 23 JA=1,NUIM 
TS(JA)= 0. 

ET(JA) = ETNEW(JA) 

2 < CONTINUE 


PRINT OUTP r ; HEADER 


IF (LONG .LT. 0.0) GO TO 44 


WRITE HEADER FOR LONG OUTPUT 


c 

WRITE (8 ,120) 

120 FORMAT (/ ’ FLUSH CODE OUTPUT’) 

CALL DATE (TAM) 

CALL TIME(TANE) 

WRITE (8, 134) TAM.TANE 
134 FORMAT (/’ A10.2X.A10) 

WRITE(8, 121) 

12 • FORMAT (/ ’ 000000000000000000000000002 '*6' ^0000000000 

< 0000000000000000000000000000000 ’ ) 

4 1 CONTINUE 


WRITE HEADER FOR TOTAL NODE CONDITION OUTPUT IF OUT ^ <' 


IF (OUT LT. 0.0) GO TO 68 
WRITE (9, 122) 

FORMAT (’ IV FI.USH1 CODE OUTPUT’) 
WRITE (9, 127) 

FORMAT (/ 

) 

WRITE (9, 124) DIAM.FLOW 
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124 FORMAT (/ ’ PIPE DIAMETER = ’ ,F6.2, IX, ’IN. ’ ,9X, ’FLOW = ’, 

♦ F7 . 2 , 2X , ’GPU’) 

WRITE (9, 125) TOTLEN.XIN 

125 FORMAT (/ ’ TOTAL PIPE LENGTB = ’ ,F8.2, IX, ’FT ’ ,4X, ’INITIAL 
. SOLID FRACTION = \F5.b) 

WRITE(9, 129) E, TO 

129 FORMAT ( / ’ PIPE ROUGHNESS = ’ ,F9. 6, 9X, ’INITIAL TEMPERATURE = ’, 
♦F7.2.1X, ’R’) 

WRITE (9, 126) 

IF (XIN .LE. 0.0) GO TO 11 
WRITE (9, 123) 

123 FORMAT(/ ’ N0DE’,3X,’ PRESSURE ’, 4X , 

♦ ’FLOW, LB/S’ ,5X, ’ FINAL SOLID * , 3X , ’ CHANGE IN X ’,4X, 

♦ ’ PRESSURE , ’) 

WRITE (9, 128) 

128 F0RMAT(9X, ’DROP, PSI’,21X, ’FRACTION’, 23X, ’ PSI’/) 

GO TO 12 

C 

11 WRITE(9 , 133) 

133 FORMAT (/ ’ NODE’ ,3X, ’ PRESSURE ’, 4X , 

+ ’FLOW, LB/S’, 5X,’ FINAL SOLID’, 3X, ’ FINAL TEMP. ,R \3X, 

♦ ’PRESSURE, ’) 

WRITE (9, 138) 

138 FORMAT (9X, ’DROP, PSI ’ ,21X, ’FRACTION ’,23X, ’ PSI’/) 

C 

12 CONTINUE 

126 FORMAT!/’ 

* * ’) 

68 CONTINUE 

C 

C REINITIALIZE VARIABLES 

C - ■--- 

C 

HP = 0.0 

WS = 0.0 

EP1(1) - PIN1 

EQUAL (1) = XIN 


BEGIN CALCULATION LOOP FOR DIFFERENT PIPE ELEMENTS 


DO 32 J=1 ,NUM 
DIAM = EDIA(J) 
WS = 0.0 
EFFP = 0.0 
LEN = ELEN(J) 

T - ET(J) 
TYPE - PTYPA(J) 


IF A PUMP IS USED, REINITIALIZE QIN AS HORSEPOWER (TYPE - PH) 
r OR PRESSURE RISE (TYPE = PP) AND CLD AS EFFICIENCY 


IF ((TYPE ,EQ. ’ PP’) -OR. (TYPE .EQ. ’ PH’)) GO TO 781 

IF ((TYPE .EQ. ’ PP ’) OR. (TYPE .EQ. ’ PH ’)) GO TO 781 
IF ((TYPE .EQ. ’PP ’) OR. (TYPE .EQ. ’PH ’)) GO TO 781 
GO TO 782 
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782 QIN EqiN(J) 

CONLD- ECLD(J) 

GO TO 783 

781 EFFP = ECLD(J) 

IF (TYPE .EQ. ' PH’) GO TO 785 
IF (TYPE EQ. ’ PH ’) GO TO 785 
IF (TYPE EQ ’PH ’) GO TO 785 
GO TO 784 


785 HP 

- EQIN(J) 

GO TO 

783 

78-1 DPP 

-EQIN(J) 

783 CONTINUE 

F 

EP1(J) 

XI 

EQUAL (J) 

7/2 

EZ2(J+1) 

71 

EZ2(.I) 

EL KM I 

NODE(J) 

JA 

.1 * 1 

SET CONSTANTS 

N 

0 

N1 

0 

N2 

0 

HFUS 

-- 25.02 

TOLD 

= ET(J) 

PI 

3 14159 

GC 

: 32.2 

DIA 

- DIAM/12 . 

ED 

E/DTA 


i ;■ 'TYPE EQ. ’ PP’) .OR. (TYPE .EQ. ’ PH’)) GOTO 786 

!F I ! TYPE .EQ. ’ PP ’) .OR. (TYPE .EQ. ’ PH ’)) GO TO 786 
IF i TYPE EQ. ’PP ’) -OR. (TYPE .EQ. ’PH ’)) GO TO 786 

.'HECK TO SEE IF HEAT LEAK INTO ELEMENT IS IN 

PIT Hit OR HTU/FT HR 

IF (HEAT .LT. 0.) GO TO 64 
GO TO 67 

64 QINTOT QIN 

GO TO 69 

i,7 QINTOT - QIN*LEN 

f,U CONTINUE 

CONTINUE 

RHOS 5 . 40 

RHOS DENSITY OF SOLID HYDROGEN, LB/CU.FT 


CHECK TO MAKE SURE THE LAST NODE WAS NOT AT THE 
SATURATION TEMPERATURE. IF SATURATION CONDITIONS EXIST 
USE THE APPROPRIATE GASPLUS ROUTINES (NORMAL LIQUID 
ROUTINES WILL NOT WORK FOR SATURATED CONDITIONS) 

iF (TSf 1) .OK. ET( J) ) GO TO 33 
GO TO 34 

KIIIII. DLPSAT (P, ID) 
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VIS1 = VLPSAT(P.ID) 

GO TO 35 

34 CONTINUE 

RHOL = DENLPC (TOLD , P , ID) 

VIS1 = VISLPC (TOLD , P , ID) 

35 CONTINUE 
C 

C USE GASPLUS TO CALCULATE LIQUID HYDROGEN PROPERTIES 
C (VISCOSITY HAS UNITS OF LB/FT-SEC) 

C 

c 

C CALCULATE FLOW PARAMETERS 

C 

C 

RHOMIX = Xl*RHOS ♦ (1. -XI). RHOL 
AREA - PI*DIA**2*0.25 
EMDOT(J)= FL0W*RH0MIX*6O . /7 . 48 
MOOT = EMDOT(l) 

VELI = MDOT/ (AREA+RHOMIX+3600 . ) 

RE = DIA+RHOMIX+VELI/VIS1 
RE1(J) = RE 

CALL THE DESIRED FRICTION FACTOR ROUTINE 

IF (XI .LE. 0.0) GO TO 17 

IP (FRICT .GB. 0.0) GO TO 17 
IP (RE .GT. 1 .08+06) GO TO 17 
GO TO 18 

18 CALL PRICT1(X1,F,RB) 

GO TO 10 

17 CALL PRICT2(F,RE,ED) 

10 CONTINUE 


SET EQUIVALENT PIPE LENGTH 


C 

C 

C IF A PUMP IS USED, THEN EFFICIENCY IS USED IN PLACE OF FLOW COEFF. 

c 

IF (TYPE .EQ. ’ PH’) GO TO 788 
IF (TYPE .EQ. ’ PH ’) GO TO 788 
IF (TYPE .EQ. ’PH ’) GO TO 788 
IF (TYPE .BQ. ’ PP’) GO TO 789 
IF (TYPE .BQ. ’ PP ’) GO TO 789 
IF (TYPE .EQ. ’PP ’) GO TO 789 
c 

IF (CONLD .GT. 0.0) GO TO 47 
GO TO 48 

47 CNLD = CONLD/F 
GO TO 49 

48 CNLD = LEN/DIA 

49 CONTINUE 

r 

CO TO 799 

r 


IF HORSEPOWER IS SPECIFIED CALCULATE SHAFT WORK, THEN PRESSURE 
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788 WS - HP*6O./(.O2356*MD0T) 

DPI A 0.0 

CNID 0.0 

QINTOT 0.0 
GO TO 790 

WS IS IN UNITS OF BTU/LB 


IK PRESSURE is SPECIFIED, THEN CALCULATE MINIMUM SHAFT WORK, WSMIN, 
ALSO KNOWN AS HYDRAULIC HORSEPOWER, THEN CALCULATE BRAKE 
HORSEPOWER , HP 


789 DP IB DPP 

DP1A 0.0 
CNLD 0.0 
QINTOT - 0.0 

WSMIN - DPP*MDOT ♦ 144 . /RHOMIX 
WSMIN IS IN UNITS OF FT-LB/HR 
HP1 WSMIN/EFFP 
WS - HP1/(MD0T*778.) 

HP HP1 ./(3600. *550.) 

GO TO 791 


!)0 


CONTINUE 

DIMA (RH0MIX*CNLD*VELI**2»F) / (2 .0*144 . *CC) 

CONTINUE 

DP IB RHOMIX* (Z2-Z1) • (32 . 2/ (GC*144 . ) ) + 

QFL0W.778 .RH0M1X/144. * EFFP*WS*778. *RH0MIX/144 

CONTINUE 

DPI DPI A * DP1B 

IF i(ABS(Z 2 -Zn LBN) .GT. 0 . 00001 ) GO TO 610 
GO TO 620 

IN HEIGHT > ELEUENT LENGTH AT ELEMENT'. 

F6 0) 

DF.LZ ABS(Z2 Zl) 

TOEMATtSx! • ’ HEIGHT L CHANOE = ' ,F7. 2, IX, 'FT. ' ,2X. 'LENGTH - ' 

F7.2, 'FT. ’) 

^FORMAT (8X/ RE-EVALUATE SELECTED FLOW SYSTEM ’/) 

CONTINUE 


VKU INITIAL FLOW VELOCITY, FT/SEC 

RF INITIAL FLUID REYNOLDS NUMBER 

VIST FLOW RATE, LB/HR 

FRICTION FACTOR 

:,n PRESSURE DROP FROM BERNOULLI'S EQUATION 


, n.T.ATK THE CHANGE IN SOLID FRACTION 

yj xi 
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30 IF (X2 .LE. 0.) GO TO 80 
C 

C CHECK TO MAKE SURE THAT THERE IS SLUSH AT THE 

C BEGINNING OF THE NODE-- IF NOT, JUMP TO LIQUID ROUTINES 

c 

P2 = P - DPI 

IF (P2 .LE. 1.3) GO TO 198 
GO TO 199 

198 WRITE(9,900) 

WRITE (9, 901) 

900 FORMAT (/ ’ EXIT PRESSURE IS LESS THAN 1.1 PSIA THE TRIPLE 
♦ POINT OF HYDROGEN’) 

90? FORMAT (/ ’ RESET INITIAL PRESSURE AND TRY AGAIN.’ ) 

GO TO 70 

199 CONTINUE 
EP1(J*1) --- P2 

IF (PI .CE. 190.) GO TO 70 
C 

C P2 = FINAL (DOWNSTREAM) PRESSURE 

C 

PAVG = (P2 ♦ P)/2.0 

C 

RHONEW = DBNLPC (TOLD , P2 , ID) 

RHOMIX = X2+RH0S * (1. -X2) * RHONEW 
VELFA = MOOT/ (RHOMIX* AREA *3600 . ) 

QFLOWA = (VELFA* *2-VELI* *2) / (2 . *GC*778 . ) 

QZ = (Z2-Zl)/778. 

HLIQI = ENTLPC (TOLD , P , ID) 

HFUS = 25.02 

HMIXI = X1*HFUS ♦ HLIQI 

HMIXF HMIXI * QINTOT/MDOT QZ - QFLOWA WS 
HLIQF ENTLPC (TOLD , P2 , ID) 

0 


<’ RHOMIX MIXTURE DENSITY, LB/CU.FT 

>' QFLOWA KINETIC ENERGY TERM, BTU/LB 

C QZ POTENTIAL ENERGY TBRM, BTU/LB 

C HLIQI = INITIAL ENTHALPY OF THE LIQUID AT T ft P , BTU/LB 

HFUS = HEAT OF FUSION, BTU/LB 
HMIXI = INITIAL ENTHALPY OF THE MIXTURE, BTU/LB 
HMIXF = FINAL MIXTURE ENTHALPY, BTU/LB 

WS = SHAFT WORK, BTU/LB 


C 

0 

(' 


I 


t 

( 

I 

( 

I 


IF (HLIQF .LE. HMIXF) GO TO 80 

CHECK TO SEE THAT THE ENTHALPY BALANCES 

QUALE =• (HLIQF-HMIXF) /HFUS 

DQUAIiP ( 144 . *DP1A) / (778 . *RHOMIX*HFUS) 

QUAL = QUALE - DQUALP 

EQUAL (J< 1 ) - QUAL 


QUALE FINAL SOLID FRACTION AFTER HEAT LEAKS HAVE BEEN 
TAKEN INTO ACCOUNT 

DQUALP SLUSH LOSS DUE TO FRICTION LOSSES 

QUAL FINAL SOLID FRACTION 
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IK (QUAL LE. 0.0) GO TO 80 

CHECK TO SEE THAT SLUSH EXISTS AT THE NODE EXIT; IF NOT, 
JUMP TO LIQUID ROUTINES 

REPEAT CALCULATIONS IN THE PROCESS OF CONVERGENCE 

RH0MIX2 QUAL*RH0S * (1. -QUAL) *RH0NEW 

RHOLI DENLPC (TOLD , P , ID) 

RHOMI QUAL*RHOS ♦ (1 . -QUAL) *RHOLI 

VIS2 VISLPC (TOLD , P2 , ID) 

VELF - UDOT/ (RH01iIX2*AREA*3600 . ) 

QFLOW - (VELF**2-VELI**2) / (2 . *GC»778 . ) 

RFY2 DIA*RH0MIX2*VELF/VIS2 

REl(.n REY2 


IF , X2 LE. 0.0) GO TO 71 

IF (FRICT .GE. 0.0) GO TO 71 
IF (RE .GT. 1 .0E+06) GO TO 71 
CAI.!, FRICT1 (QUAL , F2 , RBY2) 

GO TO 74 

CALL FR I CT2 (F2 , REY2 , ED) 
i CONTINUE 


it T i t'K EQ. ’ PH') i TO 888 
IF (TYPE .EQ. ’ PH ’) GO TO 888 
IF (TYPE EQ. 'PH ’) GO TO 888 
IF (TYPE EQ. ’ PP’) GO TO 889 
IF (TYPE .EQ. ’ PP ’) GO TO 889 
IF (TYPE .EQ. ’PP ’) GO TO 889 
CO TO 894 


*S HP.fiO /(.02356«MDOT) 

DP2A 0.0 
CO ! I "90 


US IS IN UNITS OF BTU/LB 


nr.’H 

DPP 


0.0 

*SMIN 

l)PP*MDOT* 144 . /RHOHIX 

I N 

IS IN UNITS OF FT-LB/HR 

ill’ i 

WSMIN/EFFP 

w - 

HP1 / (MD0T*778. ) 

tl'* 

IIP!/ (3000. *550.) 


T'. H>H 


i ONTINUK 


I'l’L’A 

(KH0MIX2*CNLD»VF,LF**2*F2) / (2 0*144 *GC) 

CONTINUE 


0I'2H 

RH0MIX2*QZ*778 / 144 . ♦ QFL0W*778. »RH0MIX'2, 144 


. EFFP*WS*RH0UIX2*778./144. 

CONTINUE 


DPV 

DP2A • DP2B 

ei>P2(.i* 1 ; 

DP2 

M OWI' 

(MD0T*7 . 48) /(RH0HIX2«OO. ) 
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c 


HMIXF 

QUALE 


= HMIXI ♦ QINTOT/MDOT - QZ - QFLOW - WS 
= (HLIQF- HMIXF) /HFUS 


DQUALP - (144 . *DP2A) / (778 . *RHOMIX2+HFUS) 

QUAL = QUALE-DQUALP 

EQUAL ( J » 1 ) - QUAL 
KRHO(.M) = RH0MIX2 

RECHECK TO SEE IF SLUSH EXISTS AT THE NODE EXIT 
IF (QUAL LE. 0.0) GO TO 80 

CHECK TO MAKE SURE THE DENSITY AND PRESSURE CONVERCE 

IF (ABS(DP2 DPI) .LE. 0.0001) GO TO 10 
GO TO 20 
20 DPI DP2 
N - N* 1 

IF (N GE. 400) GO TO 15 
X2 - QUAL 
GO TO 30 


10 DELRHO = ABS(RH0MIX2 - RHOMIX) 

IF ((DELRHO) .LE. 0.0001) GO TO 15 
GO TO 16 

1<> RHOMIX = RH0MIX2 

N1 N1 ♦ 1 
X2 QUAL 

IF (N1 .GE. 400) GO TO 15 
GO TO 30 


C LIQUID HYDROGEN AT ENTRANCE OR EXIT 

C ••■--- 

c 

80 CONTINUE 

<; 

(' FIRST CHECK TO SEE IF SLUSH EXISTS AT THE ENTRANCE. 

r 

IF (XIN GT. 0.0) GO TO 27 
GO TO 28 

27 MDOT - EMDOT(l) 

28 CONTINUE 

TIN - TOLD 
T TOLD 

( 

< CALCULATE THE SATURATION TEMPERATURE AT THE INITIAL NODE 

C PRESSURE. IF THE SAT. TEMP. IS EQUAL TO THE INITIAL NODE 

o TEMPERATURE ASSUME SATURATED CONDITIONS AND CALCULATE 

(’ SATURATED PROPERTIES USING APPROPRIATE GASPLUS ROUTINES 

C 

TSATl = TS(J) 

TSAT2 = TSAT(P.ID) 

IF (TS(J) GE. ET(J) ) GO TO 21 

IF (ABS(TSAT2-T) .LE. .05) GO TO 21 
GO TO 82 

(' 

■ ASSUME THE PRESSURE DROP IS EQUAL TO THAT INITIALLY CALCULATED 


i 


33 


C P INITIAL PRESSURE, PSI 

(’ PIN FINAL NODE PRESSURE, PSI 

21 DP DPI 

PIN P DP 
RHOI.IN DLPSAT(PIN.ID) 

HLIN HLPSAT(P,ID) 

RIIOL - DLPSAT(P.ID) 

GO TO 86 
v‘ CONTINUE 

r CALCULATE LIQUID PROPERTIES 

‘ (IF NO SLUSH EXISTS AT ENTRANCE OR EXIT) 

X'j o o 

QI'AI U.O 

DP DPI 

P 1 N P DP 

TSAT3 TSAT (PIN , ID) 

HLIN ENTHPC(TIN,P,ID) 

RHOLIN DENSPC(TIN,P,ID) 

RIIOL DENSPC(T,PIN, ID) 

HL1QL ENTHPC(T,PIN,ID) 

CPI CPPC(TIN,PIN,ID) 

n-J. CPPC(T,PIN, ID) 

CPI.IQ (CPl+CP2)/2. 

MOOT FL0W*RH0L*6O . /7 . 48 
WRITE(3,853) P , PIN, T.HLIQL, HLIN 
FORMAT ( / ’ P = ’ ,F5.2,2X, ’ PIN = ’,F6.2,2X,’ T 
F6.2.2X, ’ HLIQL \F7.2,2X,’ HLIN = *,F8.2) 
CONTINUE 

IF (.1 CT 1) MDOT=EMDOT(l) 

i INITIAL PRESSURE, PSI 

iP. FINAL PRESSURE 

KPM.I.n PIN 

VEl.L MDOT/ (RHOLIN*AREA*3600. ) 

V ELI N MOOT/ (RHOL*ARb,a ^600 . ) 

O'-LOAL (VEI,2**2 VELIN**2)/ (2. *GC»778. ) 

o;: (7,2 Z 1 ) / (778 . ) 

If -AT HLPSAT (PIN, ID) 

UU.IDh THE HEAT CAPACITY OF THE SOLID IN THE SLUSH 

l‘ ( V I N LK 0 0) GO TO 90 
GO I'll 01 

IlYl.l I (I O 

«"< Tii 

IF . y 1 I.L 0.0) GO TO 93 

iiMKLT (XI i *I1F1 ? S 

v,.l TO 92 
IIMKLT 0.0 

• ’ONTINUK 


■:» HUN • QINTOT/MDOT Q7, - QFLOWL HMELT WS 
ill I • M ' IIP 
:K (ilF GK HSAT) GO TO 83 
H TEMPPC(HF,PTN,ID) 


•' Hi > 3 
si, 
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IF (TL .LE. 0.) GO TO 83 
GO TO 84 

C 

C IF (T •> TSAT) ASSUME T=TSAT AND MAKE A COMMENT 

C 

83 TL - TSAT (PIN, ID) 

TS(J«1) - TL 

HF - HSAT 

HFI ( J + 1 ) = HSAT 

RHOL = DLPSAT(PIN.ID) 

84 KHONEW - RHOL 

VIS2 - VIS1 

DP2 = DP 

HLIQI = HLIN 

RLIQF HF 

REY2 - RE 

PI - P 

P2 - PIN 

F2 - F 

VELI - VELIN 

VELF -- VEL2 

({FLOW QFLOWL 

SET OUTPUT TEMPERATURE AND PRESSURE EQUAL TO INITIAL 
r TEMP, t PRESSURE OF NEXT NODE 

ET(J.l) TL 
EDP2 ( J ♦ 1 ) = DP 

IF (ET ( J + 1 ) .LT. 25.20) ET(J+1)=25.20 

r 

r CHECK TEMPERATURE CONVERGENCE 

C 

IF ((TL T) .GE. 0.01) GO TO 81 
GO TO 15 

81 N2 = N2 ♦ 1 

IF (N2 .GE. 400) GO TO 15 
T - TL 

IF (HF .GE. HSAT) GO TO 88 
GO TO 82 

H8 RHOLIN = DLPSAT(PIN.ID) 

HLIQL = HSAT 

GO TO 86 
15 CONTINUE 

IF (LONG LT. 0.0) GO TO 32 
CONTINUE 


WRITE LONG, INDIVIDUAL VARIABLE OUTPUT IF LONG > 0 

r - - ■ - 

c 

WRITE (8, 146) TYPE 

146 FORMAT (/’ ELEMENT TYPE = \A4) 

f ’ 

WRITE(8 , 100) QINTOT , FLOW , E , DIAM 
WRITE(8 ,141) OUT, LONG, XPUT.FRICT, HEAT 
WRITE (8, 142) TAMB , XRS , VAC , DIN , DOUT 
WRITE (8, 109) JA, CNLD,RHOS,RBONEW 
WRITE(8, 101) T,P,X1,LEN 
WRITE (8, 102) MD0T.RB0MIX2 
WRITE (8, 106) VIS2, HLIQI 
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WRITE (8, 103) HMIXI.HMIXF 
WRITE (8 ,107) REY2,F2,HLIQF 
WRITE!8 . 1 15) DP1A,DP1B,DP2A,DP2B 
WRITE(8, 104) P2 , DP2 , QUAL 
WRITE(8, 108) N,Nl,PINl 
WRITE (8 ,110) QUALE, DQUALP 
WRITE(8, 111) Z2,Zl ,QZ 
WRITE(8, 112) VELI,VELF,QFLOW 
WRITE(8, 113) HF,TL,HLIQL,N2 
WRITE(8 , 1 14) HSAT,TIN,HLIN 
WRITE(8 , 145) WS , BFFP , DPP , HP 

loo FORMAT (/ ’ QIN = ’ ,F10.3,1X, ’BTU/HR’ ,2X, 'FLOW ’,F9.2,1X, 

• OPM’, IX,’ E -- ’ ,F8. 6, IX, ’ DIA = ’ ,F7.3, IX, ’IN. ’) 

!11 FORMAT (/ ’ OUT - ’,F5.2,2X,’ LONG = ’,F5.2,2X,’ XPUT = ’, 

•F5.2,2X, ’FR1CT ’ ,F5.2,2X, ’HEAT ’.F5.2) 

1 42 FORMAT(/’ TAMB = ’ ,F6.2,2X, ’XRS - ’.F4.1.2X,’ VAC ’, 

•lf> 2,2X, ’DIN ’ , F6 . 2 , 2X , ’ DOUT - ’.F6.2) 

• ‘>1 FORMAT!/’ T = ’ ,F7.2,1X, ’R’,2X, * PI = ’ ,F8.3, IX, ’PSIA’ ,2X, 

•XI ’ ,F7.4,2X, ’ LENGTH = ’ ,F10.4,1X, ’FT’) 

!<>•< FORMAT!/ ’ NODE =’,I4,2X, ’ L/D = \F8.3,2X, 

• RHOS ’ ,F9.4,2X, ’RHOL = ’ .F0.4.1X, ’LB/CU.FT’) 

."2 FORMAT!/’ MDOT = ’ ,F9.2, IX, ’LB PER HR’ ,2X, ’ RHO.MIX = ’,F8.3,1X, 

• ’LB PER CU.FT. ’) 

106 FORMAT!/’ MU, MIX = ’ , 1PE14.7.1X, ’LB PER FT-S’,3X, 

• ’HI IQI - ’ , 0PF9 . 3 , IX, ’BTU PER LB’) 

i"3 FORMAT!/’ HMIX,I = ’ ,F9.3,1X, ’BTU PER LB’ ,3X, ’HMIX.F = ’,F9.3, 

• IX, ’BTU/LB’) 

107 FORMAT!/ ’ RE = ’ ,1PE14.7,2X, ’ F = ’ ,0PF10.7,2X, ’HLIQF = ’,F9.3) 
Ilf. FORMAT!/’ DP1A - ’ ,F8.4,2X, ’DP1B = ’ ,F8.4,2X, ’DP2A - \F8.4,2X, 

• ’ DP2B - ’ ,F8.4) 

» FORMAT!/’ PF - ’ ,F8.3, IX, ’PSIA’ ,2X, ’ DP2 = ’.F9.5.1X, ’PS1D’ ,2X, 
SOI, ID FRACTION ’,F10.7) 

FORMAT!, ’ N - ’ , 1 3 , 2X , ’ Nl =- ’,I3,2X,’ PIN ^ \F9 5) 

FORMAT!/ ’ QUALE ’, F10.6,2X, ’DQUALP = ’.F10.6) 

FORMAT, ’ 7.2 ’.F7.2.3X, ’7.1 ’ ,F7.2,3X, ’ QZ = ’,F9 3, IX. 

• ’DTI I LI)’) 

FORMAT!/’ VI ’ , F8 . 3 , 2X , ’ VF = ’ ,F8.3, IX, ’FT/S’ ,2X, 

’QFLOW ’ , 1PB14 .7, 2X, ’ BTU/LB’) 

FORMAT!/’ HF ’ ,F9.3,2X, ’TLIQ -• ’ ,F7.2, IX, ’R ’ , 

■JV. ’ III. IQI. ’ , F9 . 3 , 2X , ’N2 - ’,13) 

: ’ 1 FORMAT!/’ HSAT ’ ,F9.3,2X, ’TIN = ’ ,F7,2,2X, ’HUN ’.F9.3, 

IX,’ BTU/LB’) 

; If FORMAT!/’ WS - ’,P10.3,2X,’ EFPP - ’,F5.3,2X,’ DPP - \F8.3 
• ,2X, ’HP ’ ,F8 . 3) 

WRITE (8, 105) 

FORMAT (/’******♦**************#************♦****♦***«****♦*»•• 

***** ) 


: ivr hit individual node conditions or end node conditions 


"AT I N'T 
FTNKW O 
i ..I’D ( I i 0 
0 " ’0 I i , NUM • 1 
• !Di i .1 ; ,' 

IF ;FV AL(.I) LE 0.0; GO TO H5 
I MDOT • ’ V ■ 1 ) EMDOT(NIJM) 
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MPHR = MD0T/36OO. 

EPNEW = EDP2(J) + EPNEW 
DELQ = XIN - EQUAL (J) 

DELTP = PINl - EP1(J) 

IF (OUT .LT. 0.0) GO TO 36 

WR ITE ( 9 , 200) NODE ( J) , EPNEW , MPHR , EQUAL (J) , DELQ , EP 1 ( J ) 

200 FORMAT ( 13 , 3X , F 10 . 5 , 5X , F8 . 4 , 9X , F7 . 4 , 9X , F7 . 4 , 4X , F 10 . 4 ) 

IF ((RE1(J) .GT. 1 .0E+06) .AND. (FRICT .LT. 0.)) GO TO 72 
GO TO 703 

72 WRITE (9,701) 

701 FORMAT (/ ’ RE > 1.0E»06 RE LIMIT EXCEEDED ’) 

703 CONTINUE 

WRITF.(9,201) 

201 FORMAT (’ 

* - ’) 

GO TO 87 
C 
r 

85 CONTINUE 

EMD0T(NUM+1)=EMD0T(NUM) 

MPHR=MD0T/36OO . 

EPNEW = EDP2(J) + EPNEW 
IP (OUT .LT. 0.0) GO TO 36 

WRITE (9,202) NODE ( J) , EPNEW , MPHR , QUAL , BT ( J) , EP1 ( J) 

202 FORMAT (I3,3X,F10.5,5X,F8.4,9X,F7.4,6X,F10.3,6X,F10.4) 

IF (HF1(J) .GE. HSAT) GO TO 98 

GO TO 99 

98 IF (J .LT. 2) GO TO 99 
WRITE (9,204) 

204 FORMAT (’ TEMPERATURE EXCEEDS SATURATION TEMPERATURE’) 

99 WRITE (9, 203) 

203 FORMAT ( - 

’) 

87 CONTINUE 

36 CONTINUE 


IF (OUT .GE. 0.0) GO TO 56 

IF (QUAL .LE. 0.0) GO TO 76 

WRITE (9, 600) FLOW, MPHR, QUAL, DELQ, EPNEW 

600 FORMAT(1X,F10.3,5X,F8.3,9X,F7.4,11X,F7.4,7X,F10.4) 

IF ((REY2 .GT. 1.0E+06) .AND. (FRICT .LT. 0.)) GO TO 37 
GO TO 77 

37 WRITE (9,702) 

702 FORMAT(/’ RE > 1.0E+06 - RE LIMIT EXCEEDED ’) 

77 CONTINUE 

WRITE (9, 601) 

601 F0RMAT(3X, ’ - 

♦ ’) 

GO TO 56 

C 

c 

76 CONTINUE 

IF (TL .LT. 25.2) TL ^ 25.2 
WRITE(9,602) FLOW, MPHR, QUAL, TL, EPNEW 
6< •'.! FORMAT ( 1X,F10. 3,5X,F8 . 3,9X,F7 .4 ,8X ,F10. 3 ,8X ,F10. 4) 

IF (IIF GE HSAT) GO TO 78 
GO TO 79 
VK WRITE (9,604) 

6()4 FORMAT (/10X, ’ TEMPERATURE EXCEEDS SATURATION TEMPERATURE ' ) 
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79 WRITK(9,603) 

603 FORMAT (3X . ' 

’) 

56 CONTINUE 


RESET FLOW RATE AND REPEAT CALCULATIONS 


IF (FLO INC .LE. 0.0) CO TO 63 
IF (FLOW LT. FLOMAX) CO TO 60 
CO TO 70 

FLOW FLOW ^ FLOINC 
CO TO SO 
WRITE (9, 999) 

FORMAT ( ’ FLOW INCREMENT, FLOINC, WAS NOT INCLUDED 
IN THE INPUT • REVIEW INPUT DATASET’) 

CONTINUE 

END 


SUBROUTINE FRICTl (X,F,RE) 


THIS PROGRAM CALCULATES THE FRICTION FACTOR AS A 
FUNCTION OF REYNOLDS NUMBER FOR SLUSH HYDROGEN. 

THIS DATA WAS OBTAINED USING A SMOOTH, VACUUM -JACKETED 
PIPE 0.652 INCH IN DIAMETER AS DESCRIBED IN 
SINDT.LUDTKE, AND DANEY, 

MBS TECHNICAL NOTE 377, 

"SLUSH HYDROGEN FLUID CHARACTERIZATION AND 
INSTRUMENTATION," FEB. 1969. DATA WAS CURVE-FITTED 
TO OBTAIN THE EQUATIONS IN THIS PROGRAM. 

NOTES Till GEOMETRIC CURVE FIT, FG, TENDS TO GIVE 
I UK BEST FIT OVER THE ENTIRE RANGE; THIS IS EXPECTED 
A.- THE FRICTION FACTOR NORMALLY TAKES THE FORM 
I A/RE..B 

IANDARD FRICTION FACTORS GIVEN IN THE LITERATURE 
ARE. FOR TURBULENT FLOW, 

I 4*0.079*RE**- .25 AND 
F 4*0.046*RE**- .20 

THE SECOND EQUATION GIVES A BETTER FIT TO EXPERIMENTAL 
DATA. IN ADDITION, THE EQUATION OFTEN USED IS 
F.» 0.5 ^ 4.06*LOG(RE/F**0. 5) - 0.60 
THIS IS USED TO COMPARE DATA IN THE NBS DATA. 

Fi'R LAMINAR FLOW, p 64/RP 


WE FRICTION FACTOR CALCULATIONS CAN ONLY BE USED FOR 

I:! • JK06 


CF.OME PR I c M l () AX* *H) OF THE FRICTION FACTOR 
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IF (X .GE. 0.5 ) GO TO 1 
GO TO 2 

2 IF (X .GE. 0.4) GO TO 3 

GO TO 4 

4 IF (X .GE. 0.3) GO TO 5 

GO TO 6 

IF (X .GE. 0.2) GO TO 7 
GO TO 7 

I FG - 3.0775*RE**-.41267 
F = FG 

GO TO 10 

3 FG4 = 1 .94 175* RE* *-.38524 
FG5 = 3.0775*RE** .41267 

F - FG4* ( . 5-X) *10. ♦ FG5* (1 . - ( ,5-X) *10. ) 

GO TO 10 

5 FG4 = 1.94175*RE**-. 38524 
FG3 = 0.70205*RE** .309898 

F = FG3* ( .4-X) *10. ♦ FG4«(1 . - ( . 4-X) *10. ) 

GO TO 10 

7 FG2 - 0.421259*RE**-. 27205 

FG3 = 0.70205*RE**-. 309898 
IF (X LE. 0.2) GO TO 11 
GO TO 13 

II AX .2 

GO TO 12 

13 AX -- X 

12 F FG2* ( . 3-AX) * 10. ♦ FG3* (1 . -( . 3- AX) *10. ) 

10 CONTINUE 

RETURN 
END 

C 

C 

SUBROUTINE VACUUM (QRAD,TAMB,XRS,DIN,D0UT,TIN,E1 ,E2) 

REAL L,XMFP,MU 

C • .. 

0 THIS PROGRAM CALCULATES THE HEAT TRANSFER THROUGH 

C VACUUM JACKETED LINES AS A FUNCTION OF OUTSIDE 
C TEMPERATURE AND PIPE DIAMETBR FOR SLUSH HYDROGEN SYSTEMS. 

C 

c VARIABLES: 

r El EMISSIVITY OF INSIDE SURFACE, UNITLESS 

E2 EMISSIVITY OF OUTBR SURFACE .UNITLESS 

c DIN--OUTER DIAMETER OF THE INSIDE PIPE, IN 

' DOUT--INNER DIAMETER OF THE OUTER PIPE, IN 

TAMB- -AMBIENT TBMPBRATURE, R 
TIN- FLUID TEMPERATURE, R 
•' XRS — NUMBER OF RADIATION SHIELDS 

'• QRAD- -RADIATION HEAT TRANSFER, BTU/FT-HR 

* F12 -CONFIGURATION PACTOR 

( PE EMISSIVITY FACTOR 

C SIG— -STEFAN-BOLTZMANN CONSTANT 

r - 

C 

C CALCULATE HEAT LOSS DUE TO RADIATION 

C 

T2 TAMB 
T1 TIN 


CONTINUE 
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IF (XRS .GT. 0.) GO TO 15 
GO TO 16 
15 FI 2 = 10 

El IN - 1 ./El 
E2IN 1 ,/E2 
ES 0.03 

(' 

C ES ALUMINUM SHIELD EMISSTVITY 

ESIN - 1 ./ES 

CONI El IN + ESIN -1. 

C0N2 (XRS- 1 ,)*(2.0*ESIN - 1.) 

C0N3 - E2IN ♦ ESIN -1. 

FEIN - CONI ♦ C0N2 * C0N3 
FECALC- 1 . /FEIN 
FEACT = 5.0*FECALC 
FE = FEACT 
GO TO 17 

!*> FI 2 - 1.0 

A1A2 DIN/DOUT 

El IN 1 . /El 

E2IN 1./E2 

FEIN ElIN ♦ A1A2*(E2IN-1 .0) 

FE - 1.0/PEIN 
i' ;, 1 3.14159 

SIG 0.173E08 

A 1 FI. DIN/12 

TF2 T2».4. 

TF1 Tl**4. 

QRAD F12*FE*S1G*A1 * (TP2 TPl) 

QRAD HAS UNITS OF BTL/HR-FT 

20 RETURN 

c.ND 

SUBROUTINE FRICT2(F2,RE,ED) 

THIS SUBROUTINE CALCULATES FRICTION FACTORS FOR 
TURBULENT FLOW USING THE COLBROOK EQUATION. 

WRITE(14 ,201 ) RE.ED.F2 
-Ui F0RMAT(1PE15.8,2X,E15.8,2X,E15.8) 

FI 1 325'(ABS(AL0G(BD/3.7 + 5.74/RE**0.9)))**2.0 

N 0 

S'.: a 0.021 
CONTINUE 

COM (F.D/3.7 . 2.51/ (RE*SQRT(F2A) ) ) 

C0NF2 2 . *ALOG10(CON1) 

FIB 1 O/(C0NF2)**2. 

IF (ABS(F2B~F2A) LE 0.000001) GO TO 20 
F2A F2B 
N :■ N.l 

WRITE( 15 , 203) RE, N 
203 FORMAT ( 1 PEI 5. 8, 2X, 15) 

IF (N GT. 400) GO TO 23 
GO TO 30 
’<■ CONTINUE 

WRITE(9,206) RE, F2, CONI 
•-’OR FORMAT (1 PR 15 . 8 , 2X,B15. 8 ,2X , E15 . 8) 
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P2 = F2B 
23 RETURN 

END 

C 

SUBROUTINE QFLUX(QFLA.QLA) 

C 

c PROGRAM TO CALCULATE THE HEAT TRANSFER RATE 

C TO SLUSH HYDROGEN BEING TRANSFERRED IN 

C PIPING WITHOUT VACUUM- JACKETING . 

f 

C QFLUX ALLOWS UP TO 3 LAYERS OF INSULATION AROUND 

C A BARE PIPE 

<; 

C TERRY HARDY 12-16-88 

C 

DATA ID/’PH2’/ 

DIMENSION NELEM(90) , 

• TK(10, 90) , 

• INSUL(90) , HIN(90) ,H0UT(9O) ,U(90) ,QFLX(90) .AREA (90) , 

• VELI(90) ,QL(90) 


COMMON /FLUX/ EP0D(9O) ,ETI1 (90) ,ETI2(90) ,ETI3(90) , 
NNOVJ (90) ,N0JAC,NXMAT(2O,9O) , 

EXLEN (90) ,EXDIA(90) , EXT (90) ,TAMB,P,X1 ,FL0V,EPID(9O) ,K 


CHARACTER* 10 TIM , TONE , MIM , MONE 
INTEGER INSUL , NXMAT 
REAL MDOT 

C 

I - 1 
NUM = 1 
AIRV s 0 01 

c 

CALL DATE (TIM) 

CALL TIME (TONE) 

' WRITE (10, 130) TIM .TONE 

0 130 FORMAT (/ ’ A10.2X.A10) 


DO 54 I l.NUM 

IF (NXMAT (2, K) .LE. 0.0) GO TO 71 
IF (NXMAT(3,K) .LE. 0.0) GO TO 72 
IF (NXMAT(4,K) .LE. 0.0) GO TO 73 
GO TO 74 


71 

INSUL(I) 
GO TO 75 

1 

72 

INSUL(I) ■■ 
GO TO 75 

2 

7:* 

1NSUL(I) = 
GO TO 75 

3 

74 

INSUL(I) - 

- 4 

75 

CONTINUE 


54 

CONTINUE 



INITIALIZE PROPERTIES 
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r 

m 


- DO 55 I = l.NUM 

I TFL --- BXT(K) 

■ RHOL - DENLPC(TFL,P, ID) 

RHOS DNSOPC(TFL,P,ID) 

RHOMIX - Xl*RHOS ♦ (l.~Xl)*RHOL 
PI - 3.14159 
AREA (11= PI.EPID(K). *2/4.0 
MDOT - FLOV*RH0MIX*6O. /7 . 48 
VELI (I) = MOOT/ (AREA (I) *RH0UIX*36OO. ) 
CONTINUE 

DO 50 I 1 ,NUM 

DO 51 J 1 ,INSUL(1) 


MATCH MATERIAL CODE TO THERMAL CONDUCTIVITY VALUES 


IF (NXMAT(J.K) .LE. 101) GO TO 11 
IF (NXMAT(J.K) .EQ. 102) GO TO 12 
IF (NXMAT(J.K) .EQ. 103) GO TO 13 
IF (NXMAT(J,K) .EQ. 104) GO TO 14 
IF (NXMAT(J.K) .EQ. 201) GO TO 21 
IF (NXMAT(J.K) .EQ. 202) GO TO 22 
IF (NXUAT(J.K) .EQ. 203) GO TO 23 
IF (NXMAT(J,K) .EQ. 204) GO TO 24 
IF (NXMAT(J,K) .EQ. 205) GO TO 25 
IF (NXMAT(J.K) .EQ. 301) GO TO 31 
IP (NXUAT(J,K) .EQ. 302) GO TO 32 
GO TO 11 


TK(J.K) = 

0.112 

GO TO 70 


TK ; .1 , K ) 

2 717 

'III TO 70 


Tk..J,K) 

22.31 

M TO 70 


Tc. .! ,K) 

156.0 

■ . . To 70 


J.K) 

0.020 

TO 70 


TF ( J , K ) 

0 032 

i;ii TH ?u 


TK ■ .1 , K 

0.019 

Oil TO 70 


TKfJ.K) 

0.019 

Ml TO 70 


TK(.1,K) - 

0.021 

<',() TO 70 


TKlJ.K) 

0.014 

CO TO 70 


TK(J.K) 

0.015 

•'ti TO 70 


i cntinvk 


CONTINUE 


'MNTINUE 



1 'll LiiVMNC. WATER I Al CODES APPLY 
MATERIAL 
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c 
c 
c 
c 
c 
c 
(• 
( ■ 

I ■ 

r 
( ' 

(. 

c 

c 

c 

c 

c 

c 


101 

TEFLON 

102 

304 STAINLESS STEEL 

ion 

CARBON STEEL 

104 

6063 ALUMINUM 

201 

FOAM GLASS 

202 

SILICA 

203 

POLYURETHANE 

204 

POLYSTYRENE 

205 

FOAM RUBBER 

301 

FIBER GLASS 

302 

PERLITE 

THE VALUES 

OF THERMAL CONDUCTIVITY 


FROM "CRYOGENIC SYSTEMS", 2ND ED., RANDALL F. BARRON, 
OXFORD UNIVERSITY PRESS, 1985. 


CALCULATE CONDUCTION HEAT TRANSFER PARAMETERS 


DO 60 I = 1,NUM 

DR1 - (EPOD(K) - EPID(K) )/2 .0 

CONDI = (DRl/TK(l,K))/((2.0*PI*EXLEN(K)*DRl)/ 

i v ALOG(EPOD(K) /EPID(K) ) ) ) 

IF (JNSUL(I) LE. 1) GO TO 90 

R1 KPOD(K) /2.0 
R2 EP0D(K)/2.O ♦ ETIl(K) 

R3 = R2 ♦ ETI2(K) 

R4 = R3 ♦ ETI3(K) 

HI DR2 = R2 - R1 

C0ND2 = (DR2/TK(2,K) )/( (2.0*PI*EXLEN(K)*DR2)/ 
(AL0G(R2/Rl))) 

IF (INSUL(I) .LE. 2) GO TO 80 

82 DR3 = R3 - R2 

C0ND3 = (DR3/TK(3,K) )/((2.0*PI*EXLEN(K) *DR3) / 
(AL0G(R3/R2)) ) 

IF (INSUL(I) .LE. 3) GO TO 80 

83 DR4 = R4 - R3 

C0ND4 = (DR4/TK(4,K))/((2.0*PI*EXLEN(K)*DR4)/ 
(AL0G(R4/R3) ) ) 


CONTINUE 

CONI NS C0ND2 « C0ND3 ♦ C0ND4 

CONTINUE 

CONTOT CONDI ♦ CONI NS 

PAMB 14.7 

TFI, EXT(K) 

CALL SUBROUTINES FOR CALCULATION OF CONVECTIVE 
HEAT TRANSFER COEFFICIENTS 


CALL 1NHC0(X1,TFL,P, VELI(I) ,EXLEN(K) ,EPID(K) ,IIIN(I) ) 
CALL OUTHCO (TAMB , PAMB , AIRV , EXLBN (K) , EPOD (K) , HOUT (I) ) 
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HEAT I = 1 . 0/(PI*EPID(K) *EXLEN(K) *HIN(I) ) 

HEATO = 1 .0/(PI*EPOD(K) *EXLEN(K) *HOUT(I)) 

UINV = HEATI ♦ CONTOT * HEATO 

U(I) =• l.O/UINV 

QFLX(I) U(I) * (TAMB-TFL) 

QL(I) - QFLX(I)/EXLEN(K) 

QFLA -- QFLX(I) 

QLA = QL(I) 


»: QFLX - HEAT TRANSFER, BTU/HR 

< qi, HEAT TRANSFER, BTU/HR-FT 

•; HIN - INSIDE HEAT TRANSFER COEFFICIENT, BTU/HR-FT* *2- R 

<’ HOUT - OUTSIDE HEAT TRANSFER COEFFICIENT, BTU/HR-FT* *2-R 

<’ U OVERALL HEAT TRANSFER COEFFICIENT 

<’ CONDI- CONDUCTION THRU PIPE = (RO-RI)/K-A,LM 

>' COND2= CONDUCTION THRU INSULATION MATERIAL 1 

COND3r CONDUCTION THRU INSULATION MATERIAL 2 
<■ C0ND4- CONDUCTION THRU INSULATION MATERIAL 3 


WRITE OUTPUT FOR INDIVIDUAL VARIABLES 


WRITE(IO.OIO) 

■ *.!<• FORMAT ( / ’ 

« • - ’) 

< WRITE(10,91 1 ) K, EPID(K) , EPOD(K) , EXLBK(K) 

WRITE(10,920) FLOV , AIRV , TAMB , TFL 

< WRITE(10,921) RHOL.RHOS.MDOT, VELI(I) 

<• WRITB(10,912) INSUL(I), NXMAT(l.K) ,NXMAT(2,K) ,NXMAT(3,K) , 

• NXMAT(4,K) 

<■ WRITE(10,913) TK(1 ,K) ,TK(2,K) ,TK(3,K) ,TK(4,K) 

«■ WRITB(10,914) R1,R2,R3,R4 

' WRITE(10,915) CONDI, C0ND2 , C0ND3 , C0ND4 

«' WRITB(10,916) HIN(I) ,HOUT(I) , U (I) ,(JFLX(I) 

< '.-11 FORMAT (/ ’ ELEM = ’,I5,2X,’ ID - ’ ,F8.4,2X, ’OD = \F8.4, 

• 2X, ’LENGTH » ’ ,F8 3) 

’ 920 FORMAT (/ ’ FLOW = ’ ,F8.3,2X, ’VAIR = ’ ,F8.4,2X, ’TAMB - ’, 

<' • F7 . 2,2X, ’TFLUID --- ’,F7.2) 

921 FORMAT(/ ’ RHOL - ’ ,F7.4,2X,’ RHOS ^ ’ ,F7.4,2X,’ MDOT 

• F8.3.2X, ’ VFLUID = ’ ,F8.4) 

si ; FORMAT ( / ’ NINS - ’ , I4,2X, * PIPE - ’ , 1 5 , 2X , ’ MAT2 =- ’,15,2X, 

’ MAT3 - ’ ,I5,2X, ’MAT4 -=’,15) 

- Wl FORMAT!/’ K1 ’ ,F0.4,2X, ’K2 ’ .F8.5.2X, ’K3 

• FK.5.2X, ’ K4 ’ , F8 . 5) 

- 911 FORMAT!/’ R1 ’ , F8 . 4 , 2X , ’ R2 ’ ,F8 4,2X,’R3 =■. ’.F8.4.2X, 

' K4 ’ , F8 . 4) 

si , FORMAT!/’ CONI -= ’ , 1PE14 7 ,2X , ’ C0N2 ’,0PK8.5,2X, 

• ’ C0N3 ’ , F8 . 4 , 2X , ’ C0N4 \F8.4) 

<-91t» FORMAT!/’ HI ’ .F8.4.2X, ’HO ’ ,P8. 4 ,2X, ’ U \F8.5. 

<’ • 2X , QFLX ’ , F 1 2 . 4 ) 


CALCULATE HEAT TRANSFER FOR A NEW ELEMENT 


CON I NS 0.0 
CONDI 0.0 
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o o o 


C0ND2 = 0.0 
C0ND3 = 0.0 
C0ND4 = 0.0 


60 CONTINUE 
RETURN 
END 

C 

0 

SUBROUTINE INHCO (X , T , P , V , XLEN , DI A , HIN) 

(' 

i' SUBROUTINE TO CALCULATE INSIDE HEAT TRANSFER COEFFICIENTS 

I * 

CALL HYDRO (X, T,P, CP, TK, RHO, VIS) 

C 

C X 

(. T 

(' p 

C V 

C XI, EN 

C DIA 

C CP 

TK 
RHO 
VIS 
RE 

<’ PR 

C HIN 

(' 

RK DIA*RH0*V/VTS 

PR CP.VIS/TK 

IF (RE LT. 2100.) CO TO 10 
GO TO 20 

in IU NS 1 . 86* (TK/XLEN) ♦ (RE*PR*DI A/XLEN) * *0 . 333 

GO TO 30 

20 HINS .* 0.023* (TK/XLEN)* (RE* *0.8)* (PR* *0.333) 

30 HIN = HINS*3600. 

RETURN 
END 


SUBROUTINE OUTHCO (T , P , V , XLEN , DIA , HOUT) 
r SUBROUTINE TO CALCULATE OUTSIDE HEAT TRANSFER COEFFICIENTS 

CALL A IR (T , P , CP , TK , RHO , VIS , BETA) 

r 

r T AIR TEMPERATURE, R 

o P AVERAGE AIR PRESSURE, PS I 

\ AIR VELOCITY, FT/SEC 

XLEN FLOW LENCTH , FT 

DIA FLOW DIAMETER, FT 

CP SPECIFIC HEAT, BTU/LB R 

TK THERMAL CONDUCTIVITY, BTU/FT SEC R 

RHO AIR DENSITY, LB/CU.FT. 

VIS AIR VISCOSITY, LB/FT SEC 

KK REYNOLDS NUMBER 

PR PRANDTL NUMBER 

' GR GRASHOF NUMBER 

HOUT - OUTSIDE HEAT TRANSFER COEFFICIENT, BTU/HR-SQ.FT. -R 


SOLID FRACTION 
FLUID TEMPERATURE, R 
= AVERAGE FLUID PRESSURE, PSI 
= FLUID VELOCITY, FT/SBC 

- FLOW LENGTH, FT 

= FLOW DLAMBTER, FT 
= SPECIFIC HEAT, BTU/LB- R 
= THERMAL CONDUCTIVITY, BTU/FT-SEC-R 
= SLUSH DENSITY, LB/CU.FT. 

- VISCOSITY, LB/FT-SEC 

- REYNOLD’S NUMBER 

- PRANDTL NUMBER 

INSIDE HEAT TRANSFER COEFFICIENT, BTU/HR-SQ.FT. R 
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c 


C - 32.2 

RE -- DIA*RHOMIX*V/VIS 
PR = CP*VIS/TK 

GR XLEN**3*RH0**2*G*BETA* (T)/(VIS**2) 

NOTE: THE APPROXIMATION FOR BETA IS 1/TAVG = 1/T 

IF (V . I,T. 0.1) GO TO IS 
IF (RE .LT. 2100.) GO TO 10 
GO TO 20 
i ■ Hi 0.333 

PA 0.140 

HOOTS - (TK/XLEN) *CA* (GR*PR) **CB 
GO TO 30 

l*> HOOTS 1 .86* (TK/XLEN) ♦ (RE*PR*DIA/XLEN)**0.333 

GO TO 30 

J(> HOUTS - 0.023* (TK/XLEN) * (RE* *0.8) ♦ (PR* *0. 333) 

•«' HOUT HOUTS *3600. 

RETURN 

END 


SUBROUTINE AI R (T , P , CP , TK , RHO , VIS , BETA) 

DATA ON AIR OBTAINED FROM GASPLUS 

DATA ID/ ’AIR 7 
UP OPVPC(T,P, ID) 

TK TCNVPC(T, P , ID) 

HHIl DKNVPC(T,P, ID) 

VIS VISVPC(T,P, ID) 

BETA 1 0/T 

RETURN 

END 


SUHROI 7 I NE HYDRO ( X , T , P , CP , TK , RHO , VIS) 

DATA ON HYDROGEN OBTAINED FROM GASPLUS 

DATA ID/ ’PH2 V 
UP CPLPC(T,P, ID) 

TK TCNLPC(T,P, ID) 

KHOL DENLPC(T,P, ID) 

VIS VISLPC(T,P, ID) 

RIIOS DNSOPC(T,P,ID) 

RHO X.RHOS ♦ (1. X) *RHOL 


RETURN 

END 
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APPENDIX C 

NOTES FOR NASA LEWIS VAX USERS 

The following steps are required to use the FLUSH code on the NASA Lewis VAX 
system: 

(1) Create a file called FLUSH. INP as described in the users guide. 

Se/e'ai points could be of help in producing this file: 

• Namelists must always start In column 2 of your file. Therefore, the 
statements SPARAMS and $SEND must start in column 2. 

• The element variable ELE In the element description section is 14 for- 
mat Therefore, single-digit Integers need to be in column 4; two-digit varia- 
o es need to start In column 3. TYPE Is A4 format and can be placed in col- 
umns 5 to 8. Because the remaining variables (DIA, LENGTH, TEMP Q, K , and HT> 
nave FiQ.2 format, they need to be placed in lines 9 to 18, 19 to 28, et:. 

• The line with the descriptors ELE, DIA, etc., is required for the code 
to work, but because this Is A80 format, the user can vary the actual names of 
*he descriptors, if desired. 

• If the radiation heat transfer parameters are to be defined, a blank 
• ne is required between the element description section and the radiation 

neat transfer parameter section. 

2' To run the code, issue the statement 
RUN :5^lH.TLUSH]FLUSH 

>3) Output will be set to un't 9, In a file called F0R0C9.DA r , for the *•'!- 
na- outpu* conditions or the individual node output. The longer, individua' 
/arable output will be sent to unit 8 if LONG is greater than 0 

( 4 ) Note that to run FLUSH the user needs access to the GASPLUS properties 
code. To do this, simply add the following line to the file LOGIN.COM; 

ASSIGN NASAJSPT D : [ SPANK . GASPLUS ] GASPLUS . OLB LNK$LIBRARY 


I 
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TABLE 1. - ELEMENT TYPE FOR USE IN 
INPUT FILE 


Type 

Corresponding element 

PI 

Straight pipe 

BE 

Bel lows 

BA 

Bayonet fitting 

FX 

Flexible tubing 

EL 

90* Elbow 

E4 

45 # Elbow 

YV 

Y-pattern globe valve 

V A 

Standard globe valve 

AV 

Angle valve 

GV 

Gate valve 

BV 

Ball valve 

PP 

Pump, AP specified 

PH 

Pump, horsepower specified 
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I' 

I 

TABLE 2. - EXAMPLE FLUSH CODE INPUT FILE 


[Analysis of a 100-ft-1ong, 1 .5-i n .-diameter Schedule 5S Line.] 

SPAR AMS 

FL01IIN=10.0, FL0MAX=lOO.O, FL0INC=5.0, E = 0.00015, 

P = 35., XI = 0.5, 

Z1N=25. , NUM=18, 0UT=1 . , 

L0NG=-1 . , XPUT=1 . ,FRICT=+1 .0, HEAT=+1.0, 

TAMB=540 . , VAC=0 . 0 , XRS-10 . , NO JAC=0 . 

SEND 


ELE 

TYPE 

DIA, in. 

LENGTH, ft 

TEMP . , R . 

Q 

K=fL/D 

HT. ,ft 

1 

PI 

1.770 

1.000 

25.20 

0.633 

0.0 

25.0 

* ) 

VA 

1.770 

0.500 

25.20 

27 . 30 

7.58 

25.0 


PI 

1 .770 

19.00 

25 . 20 

0.633 

0.0 

25.0 

\ 

BE 

1 . 770 

1 .000 

26 . 20 

1.899 

0.0 

25.0 

r> 

BA 

1 .770 

0.790 

25.20 

14.78 

0.0 

25.0 

[) 

PI 

1.770 

19.00 

25.20 

0.633 

0.0 

25.0 

i 

BE 

1.770 

1 .000 

25.20 

1.899 

0.0 

25.0 

8 

BA 

1.770 

0.790 

25.20 

14.78 

0.0 

25.0 

y 

PI 

1 .770 

19.00 

25.20 

0.633 

0.0 

25.0 

10 

BE 

1.770 

1.000 

25.20 

1.899 

0.0 

25.0 

n 

BA 

1.770 

0.790 

25.20 

14.78 

0.0 

25.0 

12 

PI 

1.770 

19.00 

25.20 

0.633 

0.0 

25.0 

13 

BE 

1.770 

1.000 

25.20 

1.899 

0.0 

25.0 

14 

BA 

1.770 

0.790 

25.20 

14.78 

0.0 

25.0 

15 

PI 

1.770 

19.00 

25.20 

0.633 

0.0 

25.0 

16 

BE 

1.770 

1.000 

25.20 

1.899 

0.0 

25.0 

17 

VA 

1.770 

0.500 

25.20 

27.30 

7.58 

25.0 

18 

PI 

1 .770 

1.000 

25.20 

0.633 

0.0 

25.0 
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TABLE 3. - HEAT LEAK ESTIMATES 
FOR PIPES CARRYING LIQUID 
HYDROGEN 


[Estimates from CVI, Inc. ] 


Pipe 

diameter , 
in. 

Rigid 

pipe 

Flexible 

pipe 

heat leak, 
Btu/hr ft 

0.75 

0.462 

1 .293 

1 

.515 

1.545 

1.5 

.633 

1 .899 

2 

.854 

2.562 

3 

1 .080 

3.240 

4 

1.407 

4.221 

6 

1 .835 

5.505 


TABLE 4. - HEAT LEAK AND FLOW COEFFICIENT ESTIMATES 
FOR LIQUID HYDROGEN VALVES AND FITTINGS 


[Estimates from Cryogenic Energy Co.] 


Pipe 

diameter , 
1 n . 

Valve 
heat leak, 
Btu/hr ft 

Flow 

coefficient, 

Cv 

Bayonet 
f i tting 
heat leak, 
Btu/hr ft 

Y 

Globe 

0.50 

8.90 

9.0 

6.6 

7.27 

.75 

14.10 

12.0 

8.0 

7.27 

1 

21.80 

22.0 

14.0 

10.02 j 

1.5 

2 7.30 

47.0 

34.0 

14.78 1 

2 

27.30 

96.0 

50.0 

23.19 | 

i 
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TABLE 5. - EXAMPLE FLUSH CODE SHORT OUTPUT FILE 
[Node conditions not listed.] 

input Parameters 


Element 

Type 

Dia. , in . 

Length, ft 

Temp. , R 

Q, Btu/hr-ft 

K-f L/D 

i 

PI 

1.770 

1.000 

25 . 20 

0.633 

0.000 

• > 

V A 

1.770 

0.500 

25.20 

27.300 

7.580 


PI 

1.770 

IS. 000 

25.20 

0.633 

0.000 

4 

BE 

1 .770 

1.000 

25.20 

1.899 

0.000 

h 

BA 

1.770 

0.790 

25.20 

14.780 

0.000 

t\ 

PI 

1.770 

19.000 

25.20 

0.633 

0.000 

4 

BE 

1.770 

1 .000 

25.20 

1.899 

0.000 

H 

BA 

1.770 

0.790 

25.20 

14.780 

0.000 

y 

PI 

1.770 

19.000 

25.20 

0.633 

0.000 

10 

BE 

1.770 

1.000 

25.20 

1.899 

0.000 

1 1 

BA 

1.770 

0.790 

25.20 

14.780 

0.000 

12 

PI 

1.770 

19.000 

25.20 

0.633 

0.000 

i:< 

BE 

1.770 

1.000 

25.20 

1.899 

0.000 

14 

BA 

1.770 

0.790 

25.20 

14.780 

0.000 

If) 

PI 

1.770 

19.000 

25.20 

0.633 

0 . OCX) 

It, 

BE 

1.770 

1 .000 

25.20 

1.899 

0.000 

17 

V A 

1.770 

0.500 

25.20 

27.300 

7.580 

1H 

PI 

1.770 

1.000 

25.20 

0.633 

0 . 000 
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TABLE 5. - Concluded. 


FLUSH Code Output 


Pipe Diameter = 

1.77 in. 

Flow - 

10.00 gpm 


Total Pipe Length 

= 106.16 

ft Initial 

Solid Fraction - 0. 

500 

Pipe Houghness = 

0.000150 

Initial Temperature - 25. 

20 R 

Flow, gpm Flow, lb/s 

Final Solid 
Fraction 

Change in x 

Pressure 
Drop , ps 

10.000 

0.114 

0.4858 

0.0142 

0.0300 

15.000 

0.171 

0.4904 

0.0096 

0.0659 

20.000 

0.227 

0.4926 

0.0074 

0. 1 155 

25.000 

0.284 

0.4938 

0.0062 

0.1788 

20 . 000 

0.341 

0.4946 

0.0054 

0.2558 

25 . 000 

0.398 

0.4950 

0.0050 

0.3465 

40. (XX) 

0.465 

0.4952 

0.0048 

0.4509 

45.000 

0.512 

0.4953 

0.0047 

0.5690 

50.000 

0.569 

0.4952 

0.0048 

0.7007 

55.000 

0.626 

0.4951 

0.0049 

0.8462 

60.000 

0.682 

0.4949 

0.0051 

1.0053 

65 000 

0.739 

0.4946 

0.0054 

1.1782 

70.000 

0.796 

0.4942 

0.0058 

1 . 3647 

75.000 

0.853 

0.4938 

0.0062 

1 . 5649 

HO 000 

0.910 

0.4933 

0.0067 

1.7789 

H5.000 

0.967 

0.4928 

0.0072 

2.0065 

90.000 

1.024 

0.4923 

0.0077 

2.2479 

95.000 

1.081 

0.4916 

0.0084 

2.5029 

100.000 

1.137 

0.4910 

0.0090 

2.7717 
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TABLE 6. - EXAMPLE FLUSH CODE OUTPUT FOR NODE CONDITIONS LISTED 


Pipe 

Diameter - 

1.77 in. 

Flow - 10.00 

gpm 


Tota 

1 Pipe Length 

= 106.16 ft 

Initial Solid Fraction - 0.500 


Pipe 

Roughness - 

0.000150 

Initial Temperat 

ure - 25.20 R 



' - 







V ><it* 

Pressure 

Drop,psi 

Flow, lb/s 

Final Solid 
Fraction 

Change in x 

Pressure 

psi 

1 

0.00000 

0.1137 

0.5000 

0.0000 

35.0000 

■ ) 

0.00015 

0.1137 

0.4999 

0.0001 

34.9999 

.4 

0.00725 

0.1137 

0.4986 

0.0014 

34 . 9928 

l 

0.01010 

0.1137 

0.4974 

0.0026 

34 . 9899 


0.01025 

0.1137 

0.4972 

0.0028 

34.9898 

6 

0.01037 

0.1137 

0.4961 

0.0039 

34.9896 

7 

0.01323 

0.1137 

0.4949 

0.0051 

34.9868 

8 

0.01338 

0.1137 

0.4947 

0.0053 

34 . 9866 

U 

0 01350 

0.1137 

0 . 4936 

0.0064 

34 . 9865 

10 

0.01635 

0.1137 

0.4924 

0.0076 

34 . 9837 

1 1 

0.01650 

0.1137 

0.4922 

0.0078 

34 . 9835 

12 

0.01662 

0.1137 

0.4911 

0 . 0089 

34.9834 

12 

0.01948 

0.1137 

0.4899 

0.0101 

34 . 9805 

14 

0.01963 

0.1137 

0.4897 

0.0103 

34 . 9804 

If) 

0.01975 

0.1137 

0.4886 

0.0114 

34 . 9803 

111 

0.02261 

0.1137 

0.4874 

0.0126 

34.9774 

17 

0.02276 

0.1137 

0.4872 

0.0128 

34.9773 

18 

0 02986 

0.1137 

0.4858 

0 0142 

34.9701 

19 

0.03001 

0.1137 

0.4858 

0 0142 

34 . 97(H) 
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TABLE 6. - Concluded. 


FLUSH Code Output 




^AAAAAAAAAAAAAAAA 

~ ~ ~ 



aa.aa.aa.a^aaa-a 

Pipe 

Diameter - 

1.77 in. 

Flow = 15.00 

gpm 


Tola 

1 Pipe Length 

= 106.16 ft 

Initial Solid Fraction = 0. 

500 

Pipe 

Roughness = 

0.000150 

Initial Temperature - 25. 

20 R 

* ** ' - - 









N< ujn 

Pressure 

Drop,psi 

Flow, lb/s 

Final Solid 
Fraction 

Change i n x 

Pressure 

psi 

i 

0.00000 

0.1706 

0.5000 

0.0000 

35.0000 

2 

0.00032 

0.1706 

0.5000 

0.0000 

34.9997 

3 

0.01629 

0.1706 

0.4990 

0.0010 

34.9837 

4 

0.02242 

0.1706 

0.4982 

0.0018 

34.9776 

5 

0.02274 

0.1706 

0.4981 

0.0019 

34.9773 

5 

0.02300 

0. 1706 

0.4973 

0.0027 

34.9770 

7 

0.02913 

0.1706 

0.4965 

0.0035 

34.9709 

K 

0.02945 

0.1706 

0.4964 

0.0036 

34.9705 

9 

0.02970 

0.1706 

0.4957 

0.0043 

34.9703 

10 

0.03584 

0.1706 

0.4949 

0.0051 

34.9642 

1 1 

0.03616 

0.1706 

0.4947 

0.0053 

34.9638 

12 

0.03641 

0.1706 

0.4940 

0.0060 

34.9636 

13 

0.04255 

0.1706 

0.4932 

0.0068 

34.9574 

14 

0.04287 

0.1706 

0.4930 

0.0070 

34 . 957 l 

15 

0.04312 

0.1706 

0.4923 

0.0077 

34 . 9569 

16 

0.04926 

0.1706 

0.4915 

0.0085 

34 . 9507 

17 

0.04958 

0.1706 

0.4914 

0.0086 

34 . 9504 

18 

0.06556 

0.1706 

0.4904 

0.0096 

34 9344 

IP 

0.06589 

0.1706 

0.4904 

0.0096 

34.9341 
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TABLE 7. - EXAMPLE FLUSH CODE INPUT FILE FOR HEAT LEAK CALCULATIONS ON 

ELEMENTS 1 AND 9 

SPARAMS 

FL01IIN=10.0, FLOMAX= 100.0, FLOINC=5.0, E = 0 00015 
P = 35. , XI = 0.5, 

Z IN=25 . , NUM=18, 0UT=-1., 



LONG=- 1 . , 

XPUT=1 . ,FRICT 

=+1.0, HEAT 

=+1.0, 




TAMB=540 . , VAC=2 . 0 , XRS=10 . . NO JAC=0 . 




SEND 








ELE TYPE DIA, j.n . 

LENGTH, ft 

TEMP. ,R. 

Q 

K=fL/D 

HT. ,ft 

1 

11 

1.770 

1.000 

25.20 

-1.00 

0.0 

25.0 

♦ ) 

VA 

1.770 

0.500 

25.20 

27.30 

7.58 

25.0 

3 

PI 

1.770 

19.00 

25.20 

0.633 

0.0 

25.0 

4 

BE 

1.770 

1 .000 

25 . 20 

1.899 

0.0 

25.0 

5 

BA 

1.770 

0.790 

25.20 

14.78 

0.0 

25.0 

6 

PI 

1.770 

19.00 

25.20 

0.633 

0.0 

25.0 

7 

BE 

1.770 

1.000 

25.20 

1.899 

0.0 

25.0 

8 

BA 

1.770 

0.790 

25.20 

14.78 

0.0 

25.0 

9 

PI 

1.770 

19.00 

25.20 

-1.00 

0.0 

25.0 

10 

BE 

1 .770 

1.000 

25.20 

1.899 

0.0 

25.0 

11 

BA 

1.770 

0.790 

25.20 

14.78 

0.0 

25.0 

12 

PI 

1.770 

19.00 

25.20 

0.633 

0.0 

25.0 

13 

BE 

1.770 

1 .000 

25.20 

1.899 

0.0 

25.0 

14 

BA 

1.770 

0.790 

25.20 

14.78 

0.0 

25.0 

15 

PI 

1.770 

19.00 

25.20 

0.633 

0.0 

25.0 

16 

BE 

1.770 

1.000 

25.20 

1.899 

0.0 

25.0 

17 

VA 

1.770 

0.500 

25.20 

27.30 

7.58 

25.0 

18 

PI 

1.770 

1 .000 

25.20 

0.633 

0.0 

25.0 

Kadi at. 

ion 

Heat Transfer Parameters 




ELE 

0 

D. 

I . D . Ei Eo 




1 

1 . 

900 

3.334 0 . 

30 0.30 




9 

1 . 

900 

3.334 0. 

30 0.30 






TABLE 8. - SCHEOULE 5S STAINLESS 







TABLE 9. - EXAMPLE FLUSH CODE OUTPUT FOR HEAT LEAK CALCULATIONS 

ON ELEMENTS 1 AND 9 




!npu* Parameters 


K 1 em»*nt 

Type 

Dia. , in 

Length, 

ft Temp., R Q, Btu/hr-ft 

K 

f l,/l> 

) 

PI 

1.770 

1 .ooo 

25 . 20 

0.552 

0 

.ooo 

- 

VA 

1.770 

0.500 

25.20 

27.300 

-1 

/ 

. 580 

i 

PI 

1.770 

19.000 

25.20 

0.633 

0 

. (HM) 

\ 

BE 

1.770 

1.000 

25.20 

1.899 

0 

.ooo 


BA 

1 . 770 

0.790 

25.20 

14.780 

1) 

.()( )() 


PI 

1 .770 

19.000 

25 . 20 

0 . 633 

0 

. < )( H 1 

* 

BE 

I .770 

l .000 

25 . 20 

1.899 

0 

.ooo 

X 

BA 

1.770 

0.790 

25.20 

14.780 

0 

. ooo 

w 

PI 

1.770 

19.000 

25.20 

0.552 

0 

.000 

10 

BE 

1.770 

1 .000 

25 . 20 

1.899 

0 

.000 

1 1 

BA 

1 .770 

0.790 

25.20 

14.780 

0 

(XX) 

v: 

PI 

1 . 770 

19.000 

25.20 

0 . 633 

0 

ooo 

1/ 

BE 

1.770 

1 .000 

25 . 20 

1.899 

o . 

OCX) 

• i 

BA 

1 .770 

0.790 

25 . 20 

14.780 

0 . 

ooo 

; •* 

PI 

1 770 

19.000 

25 . 20 

0.633 

0 

OOO 

Ii> 

BE 

1 .770 

1 .000 

25 . 20 

1 . 899 

0 

0( )o 

! ; 

VA 

1 .770 

0.500 

25 . 20 

27 . 300 

i 

:>H0 

IK 

PI 

1 .770 

l .000 

25.20 

0.633 

0 

OOO 

Had i at i mi 

Heat 

Transfer 

Parameters 





K I 

Inner 0 . D . Outer I . D . 

E i Eo 




I 

i 

. 900 

3 . 334 

0.30000 0.30000 





i 

. 900 

3 . 334 

0 . 30000 0 . 30000 
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TABLE 9. - Concluded. 

FLUSH Code Output 


Pipe Diameter = 1.77 in. Flow = 10.00 gpm 

Tc.tal Pipe Length = 106.16 ft Initial Solid Fraction ~ 0.5(X) 


>e Roughness - 

0.000150 

Initial 

Temperature 

25.20 R 

- - 


■ - - - 



1 •'*. gpm FI 

ow, lb/s 

Final Solid 

Change in 

x Pressure 



Fraction 


Drop , }is 

10.000 

0.114 

0.4859 

0.0141 

o.ojjou 

; . ooo 

0. 171 

0.4905 

0.0095 

o.ort.v* 

20 . OOO 

0.227 

0.4927 

0.0073 

Oil 55 

25.000 

0.284 

0.4939 

0.0061 

0. 1788 

30.000 

0.341 

0.4946 

0.0054 

0 * 25S8 

35 . 000 

0.398 

0.4951 

0.0049 

0 . 3465 

40. (MX) 

0 . 455 

0.4953 

0 . 0047 

0.4609 

15 <KX) 

0.512 

0 . 4953 

0 . 0047 

0 . 5690 

50. (MX) 

0 . 569 

0.4953 

0 . 0047 

0.7007 

55. (XX) 

0 . 626 

0.4951 

0.0049 

0 . 8462 

60 . 000 

0.682 

0.4949 

0.0051 

1 0063 

65.000 

0.739 

0.4946 

0.0054 

1.1782 

70. 000 

0.796 

0.4943 

0.0057 

1.3647 

75.000 

0.853 

0.4938 

0.0062 

1 . 5649 

80.000 

0.910 

0.4934 

0.0066 

1 .7789 

85.000 

0.967 

0.4928 

0.0072 

2 0066 

90.000 

1 . 024 

0.4923 

0 0077 

2 2179 

95 (KX) 

1 081 

0.4916 

0 0084 

2 .Mil'll 

00 (KX) 

1.137 

0.4910 

0.0090 

2 VI. 
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ABLE 10. 



- EXAMPLE FLUSH CODE INPUT FOR HEAT LEAK CALCULATIONS FOR STANDARD 
INSULATION AND BARE PIPE ON ELEMENTS 3, 6, 12. 15. AND 18 


SPARAMS 

FLOMIN=50.0, FL0MAX=9O.O, FLOINLMO.O, E - 0.00015, 
P = 35. , XI = 0.5, 

ZIN=25. , NUk=18, OUT= 1 . , 

L0NG=+1 . , XPUT=1. .FRICT-+1.0, HEAT=+1.0, 

TAMB-540. , VAC=2.0,XRS=10. ,N0JAC=5 


SEND 


ELK 

TYPE DIA.in. 

LENGTH, ft TEMP. ,R 


Q 

K fL/D 

HT. , 

, ft 

1 

PI 

i 

.770 

1 .000 

25. 

20 


1 .00 

0.0 

25. 

0 

• > 

LA 

l 

.770 

0.500 

25. 

20 


27 . 30 

7 . 58 

25 

0 


PI 

i 

.770 

19.00 

25. 

■jo 


0 . 633 

0.0 

25 

0 

\ 

UK 

l 

.770 

1.000 

25. 

20 


1.899 

0.0 

25 

.0 

;> 

HA 

i 

.770 

0.790 

25. 

20 


14.78 

0.0 

25 

.0 

f> 

PI 

i 

.770 

19.00 

25. 

20 


0.633 

0.0 

25 

.0 

7 

BE 

l 

.770 

1.000 

25. 

20 


1.899 

0.0 

25 

.0 

8 

BA 

i 

.770 

0.790 

25. 

20 


14.78 

0.0 

25 

.0 

y 

PI 

l 

.770 

19.10 

25. 

20 


-1 .00 

0.0 

25 

.0 

10 

BE 

l 

.770 

1.000 

25. 

20 


1.899 

0.0 

25 

.0 

u 

BA 

l 

.770 

0.790 

25. 

20 


14.78 

0.0 

25 

.0 

12 

PI 

l 

.770 

19.00 

25. 

20 


0.633 

0.0 

25 

.0 

13 

BE 

i 

.770 

1.000 

25. 

20 


1.899 

0.0 

25 

.0 

14 

BA 

l 

.770 

0.790 

25. 

20 


14.78 

0.0 

25 

.0 

1:> 

PI 

l 

.770 

19.00 

25. 

20 


0.633 

0.0 

25 

.0 

It) 

BE 

l 

.770 

1.000 

25 

20 


1.899 

0.0 

25 

.0 

17 

VA 

l 

.770 

0.500 

25. 

,20 


27.30 

7.58 

25 

.0 

! S 

PI 

l 

770 

1 .000 

25. 

20 


0.633 

0.0 

25 

.0 

i ' * 

i ’ i ‘ ; 
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at Transfer Parameters 







ILK 


.id 


T D 

Ei 


Eo 





1 

I 

. 900 


3.334 

0.30 


0.30 





i 

. 9(X) 


3.334 

0.30 


0.30 




Miidii'-t i 

on Heat Ti 

•ansfer Parameters 






kil; 

0 

D 

T1 

T2 


T3 





3 


DIKY 

0 . 

125 0. 

125 

0 . 

125 





t) 


900 










12 


900 

0 . 

125 0. 

125 

0 

125 





13 


900 

0 . 

125 0. 

125 







18 


900 

0 . 

125 








LLK 


P 

I PE 

MAT. 1 

MAT. 

2 

MAT 3 






CODE 

CODE 

CODE 

CODE 




3 



102 

201 

20; 

3 


205 




f) 



102 









12 



101 

202 

204 


301 




1 3 



1(0 

201 

20 : 

3 






IK 



104 

101 
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TABLE 11. - MATERIAL CODES FOR 
USE IN STANDARD INSULATION 
OR BARE PIPE CALCULATIONS 


MateHal 

code 

Material 

101 

Teflon 

102 

304 Stainless steel 

103 

Carbon steel 

104 

6063 Aluminum 

201 

Foam glass 

202 

Silica 

203 

Polyurethane 

204 

Polystyrene 

205 | 

Foam rubber 

301 

Fiber glass 

302 

Perl 1 te 


61 


4 




TABLE 12. - EXAMPLE FLUSH CODE OUTPUT FOR HEAT LEAK CALCULATIONS FOR 
STANDARD INSULATION AND BARE PIPE ON ELEMENTS 3, 6, 12, 15. AND 18 






Input Parameters 


Element 

Type 

Dia., in. 

Length, ft 

Temp. , R 

Q, Btu/hr-ft 

K=f L/D 

1 

PI 

1 .770 

1.000 

25.20 

0.552 

0.000 

2 

VA 

1.770 

0.500 

25.20 

27.300 

7.580 

3 

PT 

1.770 

19.000 

25.20 

32.156 

0.000 

4 

BE 

1.770 

1.000 

25.20 

1.899 

0.000 

5 

BA 

1.770 

0.790 

25.20 

14.780 

0.000 

6 

PI 

1.770 

19.000 

25.20 

253.611 

0.000 

7 

BE 

1.770 

1.000 

25.20 

1.899 

0.000 

8 

BA 

1.770 

0.790 

25.20 

14.780 

0.000 

9 

PI 

1.770 

19.000 

25.20 

0.552 

0.000 

10 

BE 

1.770 

1.000 

25.20 

1.899 

0.000 

11 

BA 

1.770 

0.790 

25.20 

14.780 

0.000 

12 

PI 

1.770 

19.000 

25.20 

35 . 507 

0.000 

1 3 

BE 

1.770 

1 .000 

25.20 

1.899 

0.000 

14 

BA 

1.770 

0.790 

25.20 

14.780 

0.000 

ir> 

PI 

1.770 

19.000 

25.20 

37.947 

0.000 

16 

BE 

1.770 

1.000 

25.20 

1.899 

0.000 

17 

VA 

1.770 

0.500 

25.20 

27.300 

7 580 

18 

PI 

1.770 

1 .000 

25.20 

208 . 666 

0.000 

Radiation Heat Transfer Parameters 




Element 

Inner O.D. 

Outer I.D. 

Ei 

Eo 


1 


1 .900 

3.334 

0.30000 0. 

30000 


9 


1.900 

3.334 

0.30000 0. 

30000 
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TABLE 12. - Concluded. 
Bare Pipe or Standard Insulation Parameters 


Element 

Pipe O.D. 

Insulation 1 
Thickness, ft 

Insulation 2 
Thickness, ft 

Insulation 3 
Thickness, ft 

3 

1.0000 

0.1250 

0.1250 

0.1250 

6 

1.0000 

0.0000 

0.0000 

0.0000 

12 

1 .9000 

0.1250 

0.1250 

0. 1250 

15 

1.9000 

0.1250 

0.1250 

0.0000 

18 

1.9000 

0.1250 

0.0000 

0.0000 

Element 

Pipe Material 

Insulation 

1 Insulation 2 

Insulation 3 

3 

304 St Steel 

Foam Glass 

Polyurethane 

Foam Rubber 

6 

304 St Steel 




12 

Teflon 

Silica 

Polystyrene 

Fiber Glass 

15 

Carbon Steel 

Foam Glass 

Polyurethane 


18 

6063 Alum. 

Teflon 



+ + ••** + + + + ** + * + * + + + + ** + * + ** + + * + + * + + + + * + + + + + + + + + + + + + + + + + + mm 

FLUSH Code Output 




Pipe Diameter = 1.77 in 

Flow 

- 50 . 00 gpm 


Total Pipe Length = 106 

•16 ft Initial Solid Fraction 

0.500 

Pipe Roughness - 0.000150 Initial Temperature 

25 20 R 

Flow, gpm 

Flow, Ib/s 

Final Solid Change m 

Fraction 

x Pressure 

Drop , ps 

50.000 

0.569 

0.3589 

0.1411 

0.7080 

60.000 

0.682 

0.3761 

0 . 1 23ft 

1 .0145 

70.000 

0.796 

0.3884 

0.1116 

1 . 3758 

80.000 

0.910 

0.3077 

0. 1023 

1 . 7920 

00.000 

1.024 

0.4040 

0.0051 

2.2630 
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TABLE 13. - EXAMPLE FLUSH CODE INPUT WITH PUMPS AT ELEMENTS 2 AND 4 
SPAR A11S 

FL0MIN=352.3, FLOUAX-352 . 3 , FL0INC=lO.O, E = 0.00015, 

P = 25., XI = 0.5, 

Z IN-25 . , NUM=5, 0UT=1 . , 

L0NG=+1 . , XPUT=1 . ,FRICT=+1 .0, HEAT-* 1.0, 

TAMB-540. , VAC=0.0,XRS=10. ,NOJAC-0 

SEND 


KLE 

TYPE 

DIA , in . 

LENGTH, ft 

TEMP. ,R. 

Q 

K-fL/D 

HT. , ft. 

1 

PI 

1 .770 

5.000 

25.20 

0.633 

0.0 

25.0 

• ; 

PH 

1 . 770 

1 .000 

25.20 

1.690 

322 

25.0 

:i 

PI 

1 . 770 

5 . <X)0 

25 . 20 

0 633 

0.0 

25.0 

i 

PP 

1 770 

1 000 

25 JO 

2 643 

322 

25.0 

s 

PI 

1 770 

5 . (XX) 

25 . 20 

0.633 

0.0 

25 . 0 


i 


i 


U 
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TABLE 14. - EXAMPLE FLUSH CODE OUTPUT WITH PUMPS AT ELEMENTS 2 AND 4 


Input Parameters 


Element 

Type Dia. , in 

Length, ft Temp., R 

Q, Btu/hr-ft 

J* 

1 ! 

1 

PI 1.770 


5.000 

25.20 

0.633 

0.000 

O 

PH 1770 


1 .000 

25.20 

1.690 

0 . 322 


PI 1.770 


5.000 

25.20 

0 . 633 

0.000 

4 

PP 1.770 


1.000 

25.20 

2.643 

0.322 

.) 

PI 1.770 


5.000 

25.20 

0.633 

0.0(H) 

FLUSH 1 

Code Output 






Pine Diameter = 1.77 

in. 


Flow = 352.30 gpm 


Total P 

ipe Length = 

17.00 

ft 

Initial Solid 

Fraction 0.500 


Pipe Roughness ~ 0.000150 


Initial Temperature - 25.20 R 


Node 

Pressure FI 

Drop, psi 

ow, lb/ 

S 

Final Solid 
Fraction 

Change in x 

Pressure , 
psi 

i 

0.00000 

4.0056 


0.5000 

0 . 0000 

25.0000 

2 

0.78304 

4.0056 


0.4978 

0 . 0022 

24.2170 

3 

1.86134 

4.0056 


0.4893 

0.0107 

26.8613 

4 

1.07737 

4.0056 


0.4872 

0.0128 

26.0774 

5 

3.72037 

4.0056 


0.4787 

0.0213 

28 . 7204 

6 

-2.93547 

4.0056 


0.4765 

0 . 0236 

27 . 9356 
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